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AVANT-PROPOS. 


L'exposition de la méthode de Wronski pour la solution des problémes 
de mécanique céleste se trouve disséminée par fragments, dans plusieurs 
ouvrages de ce géométre, notamment dans les Prolégomènes du messianisme 
et dans la Réforme des Mathématiques (1847). Une technologie spéciale, 
l'emploi constant des considérations philosophiques, un grand nombre de 
lacunes qui laissent au lecteur le soin d'effectuer des calculs parfois trés 
compliqués, et par-dessus tout la non-énonciation du principe qui est la clef 
de tout le systéme, sont autant de causes qui, jusqu'ici, ont maintenu dans 
l'obscurité une méthode trés digne en elle-méme d'étre connue. 

M. Yvon Villarceau est le premier qui ait attiré l'attention sur la Mécanique 
céleste de Wronski (^); mais, comme is dif lui-même, il a dû se borner à 
vérifier l'exactitude des formules en partant de la méthode actuelle de la 
variation des constantes arbitraires. Il est incontestablement désavantageux 
d'appliquer des formules sans connaitre exactement la méthode suivie pour 
les établir; la connaissance de cette méthode peut conduire dans cerlains 


(*) Comptes rendus de l'Académie des sciences de Paris, 1881, t. XCII, n* 14, p. 815. Vote 


sur les méthodes de Wronski. Le mémoire, dont cette note est le résumé, n'a pas encore paru. 
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cas à des simplifications ou à des formules nouvelles, dont la légitimité ne peut 


étre établie que par la méthode elle-méme. 
Le travail actuel contient le résultat d'une étude sur l'enchainement des 


idées de Wronski, enchainement auquel nous espérons avoir donné une 
forme parfaitement claire et compréhensible. La première partie, qui forme 
le présent mémoire, renferme les relations dynamiques du probléme; la 
seconde partie contiendra l'exposé des méthodes analytiques proposées par 
Wronski pour effectuer les intégrations que contient la premiére. Ces 
méthodes, par leur nouveauté et leur généralité, ne seront pas, croyons-nous, 


la partie la moins intéressante de cet essai. 
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EXPOSITION CRITIQUE 


DE 


LA MÉTHODE DE WRONSKI 


POUR LA RÉSOLUTION DES PROBLÈMES DE MÉCANIQUE CÉLESTE. 


PREMIERE PARTIE. 


RELATIONS DYNAMIQUES. 


$ 1. Dans la méthode de Wronski, comme dans la méthode actuelle de la 
variation des constantes arbitraires, on se propose de trouver, en fonction du 
temps, les paramétres d'une conique variable sur laquelle peut étre considéré 
comme se mouvant en chaque instant, l'astre dont on étudie le mouvement, 
et l'équation de la trajectoire réelle s'obtient en remplacant, dans l'équation 
de la conique, les paramètres constants par leurs fonctions du temps. 

Telle que nous la présentons ici, elle différe de la méthode actuelle, 
principalement : l 

12 Par une plus grande généralité, c'est-à-dire que ses formules, vraies 
pour une trajectoire quelconque, constituent des relations concernant la 
dynamique générale d'un point matériel, et comprennent done la mécanique 
céleste comme cas particulier ; 

2» Par l'introduction de nouveaux paramétres variables, notamment de la 
vitesse moyenne , w, entre les vitesses extrémes sur la conique variable; du 
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paramétre p de cette conique, et, surtout, de la masse centrale sous l'action 
de laquelle elle est décrite ; 

3° Par un choix nouveau de coordonnées, — et principalement par la 
considération d'une ligne fixe dans le plan de l'orbite variable. 


$ 2. Nous rappellerons d'abord, sous forme de lemme, un théorème concer- 
nant le mouvement conique, qui sert de base à tout le développement des 
calculs, et qui nous servira ensuite à développer le principe que nous 
considérons comme la clef du systéme. 


LEMME. 


Si + G est la force accélératrice, due à la gravitation, qui agit sur une 
masse m! dans son mouvement relatif autour d'une autre masse m, c'est-à- 


dire si l'on a G — —-—", r étant la distance de m à m', on a aussi la 
relation 

t M UNE MESE cue 

où 


dt eet l'élément du temps t; 

do angle infiniment petit décrit par le rayon vecteur durant le temps dt, et 

w la moyenne entre les vitesses extrêmes de m’ à l’aphélie et au périhélie (ces mots recevant 
ici une extension qui s'entend d’elle-même) (*). 


(Nous donnons à G le signe négatif parce que cette force est dirigée de m’ 
vers m, c'est-à-dire en sens inverse du rayon vecteur de la trajectoire de m’ 
autour de m). 


(*) La relation (1) qui se déduit trés-facilement des formules du mouvement conique a une 
signification philosophique remarquable. Elle exprime que la force accélératrice G est en chaque 


instant proportionnelle à la vitesse angulaire du rayon vecteur de l'orbite. Si l'on appelle (va) la 
w(va) ` w (va)? 


T ww) or T^ 


de 
vitesse de m’ normale à ce rayon, on aura (va) = r = 元 ， 


Or, D ST est are de la force centrifuge, et le facteur 


et par Sir — G= 
t Si oú w est constant, varie pério- 
PRA A de= À =, z et u étant les vitesses extrêmes au périhélie et à l'aphélie. Ainsi, la rela- 
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PRINCIPE GÉNÉRAL (A). 


$ 3. Soit, dans le cas le plus général de la dynamique, m’ un point maté- 
riel décrivant une trajectoire sous Paction de forces motrices quelconques. 
Soit O un point donné de l'espace dont la distance à m’ soit Om! = r et 
auquel nous rapporterons, comme à un pôle, le mouvement de m’. 

Au temps (, m’ a une vitesse v dans la direction m'T de la tangente à la 
trajectoire. 

Appelons plan de l'orbite le plan variable Om'T du rayon vecteur et de la 
tangente, el décomposons les forces accélératrices qui agissent sur m’, en 
trois composantes rectangulaires, savoir : 

Dans le plan de l'orbite : , 

4° F, suivant le rayon vecteur dans le sens Om'; 

2° T, suivant la perpendiculaire au rayon vecteur dans le sens de la pro- 
jection de v sur cette perpendiculaire; 

3° P, normale au plan de l'orbite. 

Soit, en outre, o l'angle que fait, au temps £, le rayon vecteur r avec une 
ligne fixe du plan de Porbite, et supposons maintenant que l'on se donne arbi- 
trairement une fonction Ar, o, F, t) de r, 9, F, t et que w soit cette fonction. 
On pourra se représenter 10 de la maniére suivante : 

Soit w, sa valeur au temps /, pour lequel on a aussi 


TER ORT PS F =F,, 
et supposons que dans l'équation 
fr, s, F, t) =w 


on fasse 10 — w, — constante. Il en résultera pour F une fonction parti- 


tion (1) exprime une périodicité du mouvement de m' sur son rayon vecteur, c'est en quelque 
sorte un modérateur d'équilibre dynamique. Il est d'ailleurs facile de reconnaitre que la seule 
loi d'aetion en raison inverse du carré de la distance satisfait à la condition (1), car le principe 
des aires donne yat = k (constante) et par conséquent, d'aprés (1), — G =% Wronski, dans 
les Prolégoménes du messianisme, page 255, donne la loi (1) comme le résultat d'une déduction 
purement philosophique, et il y voit le principe général de toutela dynamique. Voyez dans la 


note I, à la fin de ce mémoire, ce qu'il faut penser de cette idée. 
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culiére que nous représenterons par H, et qui satisfera à l'équation 
fr, g, H, t) = u. 


Cette nouvelle équation représentera une trajectoire toute différente de la 
trajectoire considérée. 
Au temps /, on aura, dans ces deux équations, 


T = fu, 9 = 91, be w = W, PSF; : = H, 


et, par conséquent, 


Cette derniére égalité permettra de déterminer la relation qui lie les para- 
métres de la fonction H aux valeurs variables de la fonction F. On voit done 
que la fonction w est, au temps /,, le paramètre constant w, de la trajec- 
toire f(r, 9, t, H) = w,, dans laquelle les paramètres de la loi H sont déter- 
minés par la relation H, — F,. 


S 4. Application. L'application suivante rendra ce principe parfaitement 
clair. 
Soit 
dt 
f(r, y, E, t) = w = — F —. 
d; 
Il faudra déterminer quelle est la fonction H qui satisfait à l'équation 
d ` , , U 
w = — H > , Où w est supposé constant. Or, on a vu, par le lemme précédent, 
que cette relation n'existe que dans une section conique où w est la vitesse 


kw 


moyenne entre les extrêmes, et où l'on a — H = 57, E étant la constante du 
principe des aires, ou bien encore — H — Si si l'on suppose au foyer une - 
masse M — kw. 

L'égalité H, — F, deviendra A — F,, d'où l'on déduit, pour le para- 
métre M, l'expression 


M = — F. 
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Par conséquent, la fonction w — —FZ, pour une trajectoire quel- 
conque, représente en chaque instant la vitesse moyenne sur une trajectoire 
conique, tracée dans le plan de l'orbite, et au foyer O de laquelle serait placée 
une masse fictive — Fr, | 

Cette conique sera d'ailleurs supposée décrite sous l'influence de la vitesse 
réelle v au temps ( et de l'angle < de la tangente et du rayon vecteur au 
méme temps. 

Il en résulte que, si l'on parvient à calculer en fonction des trois forces 
générales F, T, P, pour un temps quelconque, la position du plan de l'orbite, 
la vitesse réelle, la direction du mouvement et le rayon vecteur, on aura 
immédiatement, d'aprés le principe précédent, en fonction des forces F, T, P, 
la détermination de tous les paramétres de la conique variable sur laquelle 
la masse m' peut étre censée se mouvoir. 

Les formules seront vraies dans le cas d'une trajectoire quelconque, mais, 
sil s'agit d'une trajectoire différant peu de la forme conique, la conique 
variable représentera d'une maniére approchée, en un temps donné, la tra- 
jectoire réelle. 


REMARQUE. 


$ 5. Le principe général que nous venons d'exposer permettrait également 
de déterminer complétement les paramétres variables d'une trajectoire donnée 
quelconque à force centrale, sur laquelle, en chaque instant, le point m’ pour- 
rait être censé se mouvoir. Supposons, par exemple, que af(») soit l'expres- 
sion de cette force centrale, a étant un paramètre. Quand a est constant, la 
trajectoire est d'une espèce particulière que nous supposons connue. Si F, est 
la valeur de la composante radiale au temps f, on aura pour le paramètre a, 
la valeur 

F, 


a = = ——— 


Ir 


r, étant le rayon vecteur au méme temps, et m’ pourra être censé se mouvoir ` 
| 2 


! rnm ven nl 
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sur la trajectoire que lui ferait décrire la force centrale 元 f(r), r étant le 
rayon vecteur de cette derniére trajectoire. Il suffit, sans doute, que la vitesse 
tangentielle soit commune aux deux trajectoires, pour que m’ puisse être 
considéré comme se mouvant à la fois sur chacune d'elles. Mais, la valeur 
partieuliére que nous venons d'attribuer au paramétre a, posséde l'avantage 
de rendre les forces centrales toujours égales dans les deux trajectoires, ce 
qui assure l'identité des trois premiers termes des développements de leurs 
rayons vecteurs en fonction du temps. C'est ce que nous expliquerons plus en 
détail au sujet des trajectoires coniques (*). 


Tel est le principe général qui nous a permis de comprendre la méthode 
de Wronski et dont nous allons maintenant développer les conséquences dans 
le cas où la trajectoire à force centrale est une conique, c'est-à-dire satisfait à 
la relation Hdt = — «d», w étant une constante et H la force radiale accélé- 
ratrice. 

Comme cette relation, facile à déduire du mouvement conique, n'est pas 
celle qu'on prend ordinairement pour point de départ, nous allons briévement 
faire voir comment, en la substituant à l'expression newtonienne G de la 
force H, on arrive plus rapidement encore à l'établissement de la trajectoire, 
et cela à l'aide d'intégrations d'une extréme simplicité. 


(*) On se demandera peut-être pourquoi nous n'avons pas simplement donné la relation H =F 
comme une condition nouvelle et arbitraire, à laquelle devrait satisfaire la trajectoire variable. 
C'est afin de rester fidéle à l'esprit de la méthode de Wronski, qui consiste à partir d'une rela- 
tion générale donnée entre la force F radiale et les coordonnées r, o, t, et à déterminer ensuite 
les variations des paramétres variables qui entrent dans cette relation, en fonction des trois 
composantes F, T, P. 
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CALCUL DE LA TRAJECTOIRE SATISFAISANT A LA RELATION Hdi = — wdọ, 
QUAND 20 EST CONSTANT. 


$ 6. Détermination de la vitesse. I est clair que la trajectoire sera tout 
entiére dans le plan de la vitesse 
et dü rayon vecteur à un moment 
donné, puisqu'on suppose ici nulles 
les forces perpendiculaires à ce plan. 

Soient (fig. 1) Ox une ligne de 
ce plan passant par le póle O et 
prise pour origine de l'angle o; 

v, el v, les vitesses de m’, paral- 


Fic. 1. 


léle et perpendiculaire à cette ligne. 
On aura : 


dv d 
CH =H cos e, <= H sin ? 


. dg 
et, puisque H = — w 77, 


D'où, en intégrant, 


vU, = — W Sin y + C, 
y. 
Da = + W COS y + Co, 


c, et c, étant deux constantes. 
La vitesse totale v, suivant la tangente, aura donc pour expression, 


w? + Ci + c? — 2w (c, sin p — c, cos y) 


et l'on voit que cette vitesse sera maximum et minimum sur deux directions 
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diamétralement opposées, correspondant aux angles o, o + > donnés 
par la relation igo, = ER Si, maintenant, l'on fait passer la ligne ori- 
gine Ox par le point du minimum de vitesse, et qu'on appelle z et u les 
valeurs absolues du maximum et du minimum de v, on aura d'abord tgo, — 0, 
ce qui donne c, = 0; faisant ensuite ọ = 0, o — dans les équations (17), 
on obtiendra, pour déterminer c, et w, les relations : 


U = W + Ca, 
Z = W — (a 


d'ou, 


Ainsi w est la vitesse moyenne entre les vitesses extrêmes de m’ sur sa 
trajectoire. 
En transportant ces valeurs dans l'expression de v ci-dessus, elle devient 


(2). . v=V wt + (w — u} — 2w (w — u) cos p =V u? + 2w (w— u) (4 — cos y) 


valeur de la vitesse réelle v en fonction de l'angle y, de la vitesse minimum u 
et dela vitesse moyenne w. | 


En fonction de ọ et des vitesses maximum et minimum, z et u, on aurait 
la formule : 


u+z\? /z—u\? (z+u)(z—u) CES T 
v 一 5 + SE AS COS 9 = == COS 9 


d'où 


= . 27 
¡ES ege LME ac Len CE Vi cos? È + z?sin? 2 
2 2° 
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VITESSES PARALLÈLE (v,) ET PERPENDICULAIRE (0,) AU RAYON VECTEUR. 


$ 7. On aura, par la décomposition des vitesses v, et v,, 


(pls  v,cosp + v, sino = — tw sin y COS y + w COS ? sin e — (w — u) sin y, 
(ve) = — v, sin a + va cos p = + w sin? y + w cos! o — (w — u) cos e 
c'est-à-dire 
(v) = — (w — uy) sin e 


(4). 


(va) = w — (w — u) cos ọ 


La vitesse (v,) suivant le rayon vecteur étant négative quand ọ varie de o 
à z, positive de x à 2z, le rayon vecteur est respectivement maximum et 
minimum pour les valeurs de la vitesse minimum et maximum. 


ÉQUATION DE LA TRAJECTOIRE. 


$ 8. En désignant par r le rayon vecteur correspondant à l'angle y, on 
aura pendant le temps di, 


dr = (v,) dt, 
rdg = (vs) dt 


d’où , à l'aide des équations (4) : 


dr. —(w—vu)singdg — — d[w — (w — v) eos el 

T w—(w—u)coss w — (w — u) cos y 
et, en intégrant, l- r = — L- [w — (w — u) eos] + constante; d’où 
(ENANA a A lef 


w — (w — u) cos y” 


en désignant par q la constante d'intégration. 


AAA FOIN Ar 
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L'équation (4!), qui donne r en fonction de w, u, o, montre que la trajec- 
loire décrite est une conique. 7 est maximum et minimum pour o = 0 et o — z. 
Soient n et m ces valeurs extrémes. On aura 


pour $220, Mit es d, 
pour 9 一 r， m(2w — u)—g. 


On en déduit 


n 9w—u n-—m w— u 


3 


m u n + Mm w 


et, comme n + m = 2a, grand axe de la conique, et que n— m = 2c, 
double de l'excentricité linéaire, on aura, pour la valeur de Pexcentricité 


numérique, 


€ 
= - = 
a 10 


En posant, de plus, p = 4 = 7 = Z, el transportant ces valeurs dans 


HA 


l'équation (4'), on obtiendra l'équation de la trajectoire, sous la forme connue 


RE. LE Qual SARA SERERE 


| A —ecose 


Le paramètre p en fonction de a, e, n et m prend les formes : 


p =n (1 —e)=m(1 =+ e), 


d’où 
=V nm (1 — e) 
nm 
qose S 
a 
et 
p= a (1 — e). 


L'introduction de la quantité e permet aussi de transformer les expres- 


(*) Le signe— au dénominateur provient de ce que nous avons pris pour origine des angles o, - 
la demi-ligne qui répond au rayon vecteur maximum. Si l'on partait du rayon minimum, on 
aurait r => Nous avons conservé les conventions de Wronski pour faciliter la vérifi- 
cation de ses formules. Il fonde le choix de l'aphélie pour origine des anomalies, sur deux 
raisons, l'une, algorithmique, qui nous a paru ne présenter aucun avantage réel (Prolég. du 
messianisme, p. 267), l'autre, philosophique, basée sur le principe de la moindre action : la 


vitesse imprimée à l'origine est moindre quand l'origine du mouvement est l'aphélie. 
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sions (4) des vitesses radiale et normale au rayon vecteur. Elles donnent 


(v) = — we sins, 


(6). 


(va) = w (1 — e coss). 
Quant à la vitesse totale v, elle devient 


Eos, E ehmilscir ni puc Mp + €? — 2e cos s. 


VÉRIFICATION DE LA LOI DES AIRES ET VALEUR DE SA CONSTANTE 
EN FONCTION DE f ET W. 


$ 9. Prenons, comme au $ 8, rdọ = (v,)dt. A l'aide des équations (6) 
cette valeur devient 
rdg = w (1 — e cos ọ) dt. 
Mais (5) donne 


. e p é s A l 
I — e cosp =, CHE 


donc 
rde = pw. dt 


et, en intégrant, f7?do = pw . t + constante. 

Mais /r2dp est le double du secteur s décrit par le rayon vecteur. On peut, 
évidemment, supposer toujours la constante nulle en faisant le secteur décrit 
nul pour / — 0, et Pon a 


IA OR QD,545170010. Z5) geli vb, 


qui exprime que les aires sont proportionnelles aux temps. 
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LOI D'ACTION DE LA FORCE CENTRALE H. 


$ 10. On a l'équation fondamentale 


de 
H = — w de 
et nous venons de trouver 
de Dr 
dt 5 
On a donc 
pw: 


CRECEN EE. iuam om ccn 


C'est la loi de Newton, quand pw? est constant. Le signe — indique que 
la force est dirigée vers le point O. Lorsque l'on considère la force cen- 
trale H comme due à Pattraction réciproque de m’ et d'une masse m placée : 
en O, on a, par le lemme énoncé plus haut, dans le mouvement relatif, 

M 
H = G = — == — z, en posant M = m + m'. On a donc, dans ce 
cas, qui est celui de la nature, pw? = m + m' = M, ce qui donne pour 
le paramétre p, l'expression 


pe = ou, réciproquement, (10) w = AE 
p 

$ 11. Les SS précédents suffisent pour démontrer que la loi Hdt = — wdọ, 
substituée à la loi newtonienne, donne, avec une extréme simplicité analy- 
tique, l'équation de la trajectoire et toutes les propriétés du mouvement, et 
qu'elle est en fait équivalente à cette loi, quand il s'agit du simple mouvement 
conique. Mais, quand il s'agit du probléme général de la mécanique céleste, 
ses avantages ne se bornent pas seulement à une simplification analytique. 
Elle établit alors, comme nous l'avons déjà fait comprendre, une relation 
directe entre la force, quelle qu'elle soit, qui agit suivant le rayon vecteur, 
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et certains éléments fondamentaux de Porbite variable, tandis que des rela- 
tions de cet ordre ne pourraient étre obtenues dans la méthode actuelle des 
constantes arbitraires, qu'aprés de longs calculs, et d'une maniére tout à fait 
indirecte. 

Nous bornerons ici cette exposition de là marche suivie par Wronski pour 
établir les formules du mouvement conique. (Voir l'indication de cette marche 
dans les Prolégoménes du messianisme, de la page 255 à la page 272.) 
Toutes les formules relatives à à ce mouvement, qu'il a déduites de celles qui 
précédent, sc sont exactes ; nous l'avons E EE elles s pourront 
facilement étre retrouvées par tout lecteur familier avec les propriétés des 


courbes du second degré. 


$ 12. Il nous sera facile, maintenant, gráce à notre principe général, de 
passer du cas de w constant dans la relation Fdt = — wdy, au cas de w 
variable sous l'influence des forces quelconques F, T, P. Mais, comme le plan 
de l'orbite et les lignes-origines dans le plan de l'orbite, deviendront générale- 
ment variables, il convient de s'entendre auparavant sur la maniére dont 
seront définies les coordonnées de la conique idéale sur laquelle m’ sera censé 
se mouvoir en chaque instant. 

Par le póle arbitraire O (fig. 2) ou, ce qui revient exactement au méme, 
par le centre m de la masse fictive pla- 
cée au foyer O de la conique idéale, 
faisons passer un plan fixe, et dans ce 
plan, par ce méme centre, une ligne fixe, 
mA. Le plan variable de Porbite de m’ 
coupe le plan fixe suivant la ligne des 
nœuds m, en considérant ici le côté 
du nœud ascendant, et la position de 
cette ligne des nœuds est définie par 
l'angle AmQ = y, longitude du nœud 
ascendant. Le plan mtm de l'orbite fait 
avec le plan fixe un angle », qui est son inclinaison, et quand y et » sont 


connus, la position du plan de l'orbite est entiérement déterminée. 
3 


Fic. LS 


r 


omnt —M 
à 
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Soit maintenant, dans le plan de l'orbite, une ligne fixe C, pour servir 
d’origine aux angles mesurés dans ce plan. La position de cette ligne sera 
déterminée par l'angle y = CMQ + y. 

Enfin, si 9 et « sont les angles que font avec cette ligne fixe, le rayon 
vecteur mm' — r, et le demi-grand axe mB pris du côté de l'aphélie, on 
aura la relation 


(14) e E D B . D D E . . . D ?一 和 一 0， 


en appelant, comme ci-dessus, y, l'angle du rayon vecteur avec le demi-grand 
axe du cóté de l'aphélie. 

y et » déterminent le plan de l'orbite; 

y et «, la position du grand axe dans ce plan. 

En y joignant ce grand axe lui-méme, 2a, et Pexcentricité e, on aura 
toutes les données nécessaires pour fixer les dimensions de l'orbite et son 
orientation dans l'espace. 


8 13. Abordons maintenant le probléme général que nous nous sommes 
posé. 


PRINCIPE (B). 


Dans les calculs qui suivent nous aurons plusieurs fois à évaluer la diffé- 
rentielle premiére du rayon vecteur par rapport au temps, due à l'action des 
forces accélératrices ; il importe de rappeler, pour éviter tout malentendu, le 
principe de dynamique trés connu en vertu duquel cette différentielle est 
toujours nulle. En effet si, au temps t, la vitesse du mobile sur la trajectoire 
est égale à v, le chemin parcouru dans la direction du rayon vecteur pendant 
le temps dt, sera dr = v coss - dt, et ce chemin ne dépendra des forces aecé- 
lératrices que par la vitesse v et l'angle «. Ainsi, par exemple, si pour plusieurs 
mobiles m,, m,... mus par des forces différentes et décrivant des trajectoires 
différentes, on a en un méme moment, ou en des moments différents, 


Uy = V = "+ —=V, T = Ha MAKE 
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Vis Vars, 51) Baz ee Étant les vitesses réelles et les angles de leurs directions et 
de leurs rayons vecteurs, il y aura aussi égalité entre les variations dr,, dr, ... 
de leurs rayons vecteurs pendant le temps dí. Les différentielles d'ordres 
supérieurs seront seules différentes en général. De ce principe, nous tirons 
les remarques suivantes. 


Si m,, Ma, ... décrivent des coniques autour d'une masse M, avec des 
vilesses v et des angles + égaux, quels que soient d'ailleurs leurs rayons 
vecteurs, les variations dr,, dr,,... de ces rayons vecteurs aux temps corres- 


pondants seront indépendantes de la masse M, ou des forces attractives 
M M 


esse ses 
AEN 


11. 


Si un mobile m’ décrit une trajectoire sous l'action de certaines forces, et 

qu'à un moment donné, on lui applique une nouvelle force motrice finie, — 
si l'on veut, une force perturbatrice, -— la variation du rayon vecteur pen- 
dant l'élément de temps n'en sera pas altérée, et, par conséquent, la variation 
due à l'influence de la foree perturbatrice sera nulle. 


Corollaire. .La variation du rayon vecteur d'une trajectoire quelconque est 
égale à celle du rayon vecteur de la conique idéale variable sur laquelle le 
mobile est censé se mouvoir. 

À ces théorémes nous joindrons encore le suivant : 


Les paramétres d'une conique sont complétement déterminés par la 
vitesse v, l'angle +, le rayon vecteur r et la masse attractive M. Dans la 
conique variable, nous avons vu que cette masse attractive, qui n'est alors 
que fictive, est égale à — F7?, F étant la composante radiale des forces accé- 
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lératrices. Il en résulte que , si — Fr? ne changeait pas, les paramètres de la 
conique resteraient constants (^). Les variations de v, +, r dues à la force F, 
pendant le temps dt, sont celles que subissent ces variables dans la conique 
méme, puisque l'action centrale est égale à Li — F et, par conséquent, 
les variations correspondantes des paramétres de cette conique sont nulles. 
Les variations réelles de ces paramétres se réduisent donc à celles qui pro- 
viennent des variations de la masse fictive M, c'est-à-dire de F7?. Par consé- 
quent, la force F n'agit sur les paramètres de la conique qu'en faisant varier 
la masse fictive M, et, si l'on cherche les variations de ces paramétres , il faut 
y considérer comme nulles les variations de v, « et v en tant quelles pro- 
viennent de F. | 


8 14. Ceci posé, venons à l'établissement des formules du mouvement dans 
la trajectoire, sous l'influence des trois composantes quelconques F, T, P, et 
occupons-nous d'abord des deux premiéres F, T, qui agissent dans le plan 
variable de l'orbite. 

Cherchons en premier lieu les expressions des deux paramétres variables p 
et w, qui sont les facteurs principaux dans l'expression de la masse centrale 
fictive M et dans le principe des aires. 

Nous avons trouvé pour w — constante , dans le mouvement conique, 


(I0): ie A UE ROO AC REED RN L 


Nous aurons donc ici (principe général A), quand p et w seront variables, 
LEE aru do TEMP E E 


Quant à l'expression générale de la constante du principe des aires, 
pw (éq. 8), il est clair que la force centrale F n'exerce aucune action sur sa 
variation. Mais il en est tout autrement de la force T, qui fait varier directe- 
ment la vitesse (v,) = w (1 — e cos?) =+? (éq. 5 et 6), perpendiculaire au 


7 
rayon vecteur 7. 


(*) Il n'est évidemment pas question ici des forces T et P. ` 
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La variation d(v,) de (v,) pendant le temps dt, ne sera plus seulement due 


à la variation de 7 dans l'orbite, mais aussi aux variations des paramètres w 7 


et p pendant ce méme temps. On aura donc, en représentant par d(v,) f 
la variation de (v,) sur l'orbite conique, pour la variation totale de cette d 
vitesse, 
d. wp 1 
(ve) = d(v;) + Tdt = ——— + wpd—. d 
GEN ui 4 i Q e 


Mais, comme nous l'avons déjà démontré (principe B), la variation de r 
provenant des forces accélératrices est nulle, ou, si l'on veut, un infiniment 
petit du second ordre; par is un dr = dr, variation de r pn l'orbite 
conique. On en déduit : d(v,) = wpò = et, par conséquent, Tdt = ; d'où, 
en intégrant, 


n8) vd) 424508 di 5 op e bie + d Trdt, 


(wp), étant la valeur de wp au temps (t). 

Cette équation, en prouvant que, pour T = o, wp est constant, vérifie le 
principe des aires pour le cas d'une force centrale F quelconque. 

La combinaison des équations (12) et (13) donne alors, pour la détermi- 
nation des deux variables fondamentales w et p, les expressions : 


E Epa 
E EE 
(wp)o rd Trdt 
D 
| wm e. wu | 
ENS CLAUSA uox o odere ese. 


— Fr* 


Quant à la vitesse perpendiculaire au rayon vecteur, on aura, au temps f, 
par le principe des aires, 
pw == (vs) r. 
D'où : 
t 
(wp) +, "i Trdt 


r 


(4 6) COR S em Dos T (va) 一 一 


Wi 
(29 


E 
7 


+ NP d. 
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On en déduit également (v,)dt = dr tgs, = étant l'angle de la tangente et 
du rayon vecteur; d'oü : 


t 
(wp)o n T Trdt di 


Lt UN H DEN e en € 


et, pour la vitesse réelle v, sur la vraie trajectoire, 


KA ME CAES RUPEE TENIS Em 


sins étant donné par l'équation (17). 

Ainsi, si l'on connaissait le rayon vecteur » et la force T en fonction du 
temps, les équations (17) et (18) donneraient l'angle « et la vitesse v, c'est-à- 
dire les deux éléments qui, joints à ce rayon 7 et à la masse fictive — Fr?, 
déterminent complétement, pour le temps /, les dimensions et l'orientation de 
la conique variable. | 

Langle de position 9, compté à partir de la ligne fixe du plan de l'orbite, 


et qui détermine done, à chaque instant, la position absolue du rayon 
vecteur dans ce plan, se déduit également sans peine des lois fondamentales 


(4) . . H B . . 2 . e e — Fdt — wd», 
EV oM M teg — Fr! 一 pu; 


en remarquant que, dans (1), de est le déplacement angulaire absolu du 
rayon vecteur et que, par conséquent, 


d; = d 
On trouve aisément l'expression z 
A t w pu fe 
GE 5 一 Oo+ f Da (Le 
() 
APE De gl TN EL E 


où (6), est l'angle du rayon vecteur et de la ligne fixe de l'orbite au 
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temps (1). Dans cette expression, pw peut étre remplacé par l'expression (13); 
ce qui donnera 


' dt ' dt j 
(20) — .. . E (pto)g Sé ^ 天 * f- Trdt. 
( (0 


IU) 


$ 15. Actuellement, puisque l'on possède les expressions de la vitesse 
réelle v, de l'angle < et de la masse attirante fictive — Fr? placée en m, en 
fonction des forces accélératrices et du rayon vecteur r, on en déduira trés 
facilement, par les formules du mouvement conique, les expressions du demi- 
grand axe a, de l'excentricité e et de la longitude « de l'aphélie de la conique 
variable. 


Grand axe 2a. Les équations du mouvement conique donnent facilement 
la relation : 


TENIS T) 

n) — 9pu* — rv?! 
pu? et v sont donnés par (12) et (18). 
En fonction de F, r, v, on aura donc 


2 
a. m. o 
2Fr + v* 
On peut aussi donner á a une autre forme; des équations (13) et (18), 
on déduit 


pw’ 
Es Ee 
d’où 
r $ 
Ee cs oes usb a er D ah 
BEE 
r sinto 


Cette dernière équation donnera le demi-grand axe quand, le paramètre p 
étant connu par l'équation (15), on aura aussi r et «. En remplaçant d'ailleurs 


24 EXPOSITION CRITIQUE 


dans l'équation (21), pw? et v par leurs valeurs générales (12) et (18), on 
obtiendrait finalement l'expression générale 


: Fr? 
(39) . . , . T me AS 
dr\2 | co +S Trat | 
2Fr + (5. (CASUM Ip ns 
Excentricité e. La relation p = a(4 — e?) donne immédiatement 


(23) e = V 1 —£, où tout est connu en vertu de (15) et de (22). 
On obtiendrait e en fonction de p, r et e par la transformation suivante : 


en gf (B) 1 
r? sin?s r r/ sina 
D'où 
"it 
OE gu NO TEE E 
r r) sin?c 


Longitude a de l'aphélie. L'angle D — « que fait le rayon vecteur avec la 
ligne des apsides, est donné dans le mouvement conique par l'expression 


P 
r 


p 
t AE 
e 


et, en y remplacant p par (15) et tgz par (17), 


y dr 
| + Z trar |e 
GN) a TES De EE I 


多才 一 q 一 二 一 一 一 一 一 一 


KZ + F Tu | + Pn 
N 


(25) Le" e TA EN LA 5 O tg(9 — a) = 


$ 16. Les équations précédentes contenant des intégrales qui supposent 
connue la valeur de », en fonction du temps, il faut encore connaitre cette 
fonction ou, si l'on veut, l'équation dela trajectoire. — 

Or, il est clair que r est, au temps £, le rayon vecteur d'une conique dont 
les éléments sont donnés par les équations précédentes, et qui correspond à 
l'angle 9, également donné par l'une de ces équations. 
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On a done, pour l'équation de la trajectoire dans son plan, 
EEUU v you GR f [E —-e0$ (0 — all =p, 


et cette équation, en y remplaçant e, 9, «, p par les fonctions calculées (25), 
(49), (25), (15) des composantes F, T, de r et de t, représentera implicite- 
ment toutes les trajectoires que peut décrire dans un plan un point matériel. 

Nous indiquerons plus loin la maniére de se servir de ces relations pour en. 
dégager les inconnues dans le cas de la mécanique céleste (*). 


$ 17. Avant de passer à l'étude des effets de la composante P, normale au 
plan de l'orbite variable, il nous reste à exposer d'autres relations, existant 
entre les paramétres variables, et qui peuvent également servir à les déter- 
miner en fonction les uns des autres. 


Excentricité e. Reprenons l'expression (24) et faisons y varier les variables, 
en remarquant que, d’après le principe (B) et le théorème HI, dr est nulle, 
et qu'il faut également regarder comme nulles les variations E v et < qui 
dépendent de F (**). 

On obtient ainsi : 


A iti 1 j T À 
D ER. 
i sino Ar r r) sis — 


(*) L'angle q est donné par son expression (19) en fonction de r et de t, mais d'une manière 
implicite. Les équations (19) et (26) constituent-done deux relations entre r, ® et L mais ne 
résolvent pas mieux, en réalité, le problé trois corps que les méth ¿ja connues. La 
prétention de Wronski (Réforme des mathématiques, p. c.i) de donner la solution absolue du 
probléme des trois corps est done mal fondée, Si Pon transporte, dans l'équation conique, les 
expressions des paramétres variables trouvées par la méthode actuelle des constantes arbitraires, 
on arrive à une équation de la trajectoire réelle, aussi bien gue dans la méthode de Wronski. 

(**) Notre théorème III montre Pinexactitude partielle de l'expression dont Wronëki se sert, 
pp. CXLVII et CXLIX, quand il dit que la composante radiale n'exerce aucune influence sur les 
variations de e et a, Jl faut remarquer, en outre, que la composante radiale dont il parle, n'est 
pas notre composante E, mais seulement ce qu'on entend aujourd'hui par la composante per dl 
batrice suivant le rayon vecteur, 11 se place, dés l'abord, dans le cas de la mécanique céleste. 


4 
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où l’on a, en vertu des principes rappelés, 


T cos% étant la composante de T perpendiculaire à la tangente, et qui agit 
donc pour faire tourner cette tangente de l'angle dz pendant dt. (Voyez plus 


loin pour la démonstration trés élémentaire de cette expression le $ 33.) 
On a , d'ailleurs, 


Trdt = d(pw) (éq. 15) cet rvsins — pw. 
D'où : 


Ce d(pw) 
pw 


Sins Coso. 


Par conséquent, (27) deviendra 


— à mg — — — 


PAK) o 
La LE dp a 4 p dp (2) coss d(pw) 
, w 


- —— . —— sins COST 
1 sino rv tr r} sins p 
? d 2 
MESES, uc LR 
row r' sine p r 


pv. 
gr Lr EC Loa: 161 bn) 


Or, on a, par le mouvement conique, 


i rdo 1 — ecosọ 
-—— esing ? 
rU — e cos ọ}) =p 
D'oú | 
D 
2 T 


an 
e sin! + D 
» 


1 2 
(2) Ka (P — sine) — et sint + (2) — P — eint e PP Li) mess — P ens. 
r/ sin*s Vr r T rr r 
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Ainsi, (28) devient 
ede — G sin*p —te cos ) L e sin°p E 


l d 
de = E due — cosp) m e sin*e d(pw) 


et, en intégrant, de (1) à £ et appelant (e), la valeur de e pour £, 


ou , comme 


La 
(4) 


t ; A / 
A Pn tum er fosa + eint (o — E K - 169 kad 


$ 18. Longitude « de l'aphélie. On a la relation 


p 
— = | — e cos (D — a). 
$ e cos( a) 


Nous rappellerons ici, comme au sujet de l'excentricité, le principe (B) et 
le théoréme III. Nous remarquerons, en outre, que si l'on fait varier toutes 
les quantités de l'équation précédente, on peut faire abstraction de la partie 
de ces variations qui auraient lieu dans l'orbite supposée invariable, attendu 
que ces parties des variations donnent, dans les deux membres de l'équation, 
des valeurs égales. 

Comme ces variations dans l'orbite proviennent de la variation de l'angle 9, : 
cela revient à considérer cet angle ® comme constant. 

En vertu de ces principes, nous aurons d? = o, dr — 0 et, par conséquent, 


d 
D a esin (8 — a) de + cos(® — a) de = o 


28 . =- EXPOSITION CRITIQUE ` ` 


et, puisque 


eos (6 — a) = 


Ge e sin (D — a) da 一 


fs me 167 


komm sooo 
` cesin (6 — a) e 


D'où y 


et, à cause de l'équation (29), employant pour simplifier o = ®—«, ‘+ 


D à P 

1 ro ‘p Y p dp 

de = —— {© .= 0089 —— esintp — — — - — 

e sin ($ — ol e r p € 10 rpi 

P o 
s IU re pear ae: DEL 
esin(b — a) r/ C p p 10 
r 
Or, 
p 
1 Le LONA 
donc 
p 
CS r ert sea B. 
e sin ($ — a) p B sw 
T 
ou, comme 
q A dw —d(pw) dp 
w pw p’ 
p- | 
T 


“— esin (9 — a) j p p pw 


pes E + | 


— sin? t e à R E E . 


U 
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Tenant compté enfin de la relation © + ecoso = 1, remplacant- + par 
® — a, simplifiant et intégrant, il viendra ` 


ABUS 6 ves «— (Que f sino — dlar — a gpw) y 2 


(2) 


Les équations (29) et (30) donneront respectivement e en fonction. de a, 
et « en fonction de e, quand on aura déjà les expressions des paramétres fon- 
damentaux p et w (car est également fonction de p et w par l'équation (19) ). 
^ Afin de ne pas compliquer cette exposition, nous renverrons à la note II 
quelques considérations sur une formule exponentielle, présentée par Wronski 
pour représenter la trajectoire dans son plan et que nous ne croyons pas 
exacte. 


$ 19. Les formules précédentes contenant toutes les conséquences de 
l'influence des forces F et T situées dans le plan de l'orbite, il ne nous reste 
plus qu'à évaluer celle de la force P normale à ce plan. Cette derniére force 
agit, à chaque instant, pour faire varier l'orientation du plan de Porbite. Nous 
allons, d'abord, déterminer l'angle do dont elle fait tourner ce plan autour du 
rayon vecteur 7 pendant le temps dí, et de la valeur de dp nous déduirons 
ensuite, trigonométriquement, les variations dy, dy, des angles qui fixent la 
position de ce plan. (Voyez § 12.) 

Soient Av (fig. 3) la direction de la vitesse réelle v, r le rayon vecteur, < 
l'angle de v et r. AB étant l'élément vdt 
décrit pendant le temps dí et BC, une per- 
pendiculaire à la direction du rayon vecteur, 
on aura CB = vdt sin w. | | 

La force P fera au-bout du temps di 
décrire à B un angle de autour de C dans le 
plan (CB, P), et cet angle sera celui dont le. plan (r, v) a tourné dans le 
méme temps autour de ». Or, pendant le temps di, le point d'application B 
de P a subi par rapport à A un déplacement relatif P(dt)? (c'est-à-dire 


Fic. 3. 
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dv'dt = Pdt : dt; en effet, le point A et le point B sont animés tous les deux 
au temps / de la méme vitesse acquise v' sons l'influence de P). On a 
donc 


P(dt? = dọ - CB = de vdl sine, 


Si, maintenant, e, p, et « sont les paramétres de la conique variable au 
temps 7, on aura, en posant 6 — a = 5, 


[or = n GN ar l'équation r Lal A RES 
dr  resin(® — a) P 1 4 — e cos(® — a) 


D'où l'on déduit 


dy 2 “y? 
1 +») — + €? sin? —_—— 
A \dr V SMET hity S < cos (9 — x) 
sino ; dy NT p Se À — e eos( —.«) 
dr 5 


| 
I 


Pdt VA + e — 2e cos(b— a) , e: H 
ME EE desen at TI We 
Me 


Au lieu d'introduire l'excentricité e, et les angles ® et a, on pourrait 
également se servir de la formule qui donne plus directement tgz (éq. 17). 


"n 
el, par conséquent, LP de Y 


(*) Wronski a écrit dp — = Le facteur par lequel il faut multiplier cette expression 
inexacte est celui que M. Y Villarceau a également donné dans les Comptes rendus de PAca- 
démie des scienees de Paris, dans sa Note citée sur les méthodes de Wronski. 
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$ 20. L'expression de dp étant obtenue, il est facile d'en déduire les 
expressions de » et y. Reportons-nous (fig. 4) 
à la figure du $ 12 et supposons que le plan 
de l'orbite ait tourné, autour du rayon vec- 
teur mm!, d'un angle p. L'intersection de 
ce plan avec le plan fixe deviendra m£, n 
deviendra ;/ et, Ai, 

Abaissons Parc de grand cercle m!! per- 
pendiculaire à AQS’. Nous aurons. 


Fic. 4. 


tg m'l — sinl - tg y, 


d’où 


nlo 
cos 152 tgy'd IQ + uio =:0, 
cos? y - 


On a aussi dans le triangle £262^m/ : 


sin RQQR  sinm'S2 


sin p siny’ ` 
Pour » = zl, 662 = — d - 182, p = dp, ces relations deviennent : 
in/ 
cos (£2 - tg» d - 152 +- d dy — 0, 
COS Y 
sin m'Q2 


dy — — p0LIK2 siny cosyd- 162, — d- IQ = dp. 


sins 
et, par suite, 
dy = cot 152 sin m'£2 cosy - dp. 


Mais 
tg £2 = tgm'£2 - cosy, 
donc 
(33) AS a VES NS LT e o qi). de. 
Y 


tg m'6 - cos 
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On aura, par les mémes formules, 


sin m'Q2. 


(54) ISO GHORBOGUGUS T3 dy = $262 = + de. 


sin y 


Langle oi. d’après les conventions adoptées au $ 12, est égal à 
x —4 +0, On aura done, en intégrant ES et (34) du temps (1) au 
temps /, 


tM MAC oc (do TJ cos — y + 9) de, 


| | P M 
COMER MORE en f RL A TA 


sin y 


Quant à l'angle y, qui détermine la position de la ligne fixe du plan de 
l'orbite, on aura pour le déterminer, la relation : 


A 一 由 +0 = HUK. 
Dans la rotation, l'angle 9 reste constant; done .. : 


dy = dy + d- mQ.. 


Or, v ee 
sin m'l = sin m'Q2 - sin y. 
Don: ` oW | 
d. MQR = — 1gm'$2 cot; : dy, 
— — tgm/'S coty cos m'S2 - de, 
= — sin m'62 cot : dp... 
Done 
1 A — cosy -. Gas pum 
dx = sin m'£2 - e 一 eu de = O: “de = sin m'Q2 tg2 : do 
sin y sin y 2 
et, en ee 
t ` 
: 4 = 
Bases Y PAL de era ef sin (x — y +0) tg% de. 


.的 
Dans toutes ces expressions (33), (36), (31) (*), dea la valeur (31) ou (38). 


(*) Données sans démonstration par Wronski, p. 291 des Prolégoménes du messianisme. 
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EMPLOI DES FORMULES GÉNÉRALES. 


$ 21. Le résultat immédiat des relations dynamiques établies dans les 
paragraphes précédents, consiste dans la détermination des paramétres 
variables d'une trajectoire arbitraire à force centrale, sur laquelle un point 
matériel décrivant une trajectoire donnée quelconque, pourrait être consi- 
déré comme se mouvant en chaque instant. En effet, le principe général (A) 
fournira toujours pour le paramètre de cette force H elle-même, une valeur 
(ou, si cette force H renfermait plusieurs paramètres, une condition entre 
ces paramètres) donnée par l'équation générale H — F. La force H étant ainsi 
connue et les expressions générales (17) et (18) de v et < subsistant toujours, 
on n'aurait qu'à exprimer les paramétres 4, l, m... de la courbe variable 
choisie, en fonction de v, +, H et r, ce qui sera donné immédiatement par la 
nature connue de cette courbe. 
= C'est ce principe général que nous avons appliqué plus haut, en choisis- 
sant pour trajectoire variable une conique, les intégrales contenues dans les 
expressions des paramétres variables étant calculées au moyen des coordon- 
nées de la trajectoire réelle. 

Dans la question présente, l'emploi le plus intéressant de ces relations 
dynamiques consiste à déterminer d'une maniére de plus en plus approchée 
la trajectoire d'un point matériel, quand on connait déjà cette trajectoire 
par une premiére approximation. Il suffit pour cela de prendre pour trajec- 
toire variable la courbe méme donnée par cette premiére approximation, et 
d'en calculer les paramétres variables par la méthode générale exposée plus 
haut pour une trajectoire variable quelconque. La trajectoire variable ainsi 
choisie fournira pour toute époque déterminée /, une représentation géomé- 
trique plus exacte de la trajectoire réelle. 

D'aprés cette régle, quand il s'agit des mouvements célestes, dont la pre- 
miére approximation est une conique, la trajectoire variable conique s'impose 
done d'elle-méme. On exprimera d'une maniére approchée, en fonction du 
temps, les forces F, T, P et le rayon vecteur et l'on pourra alors effectuer 
les intégrations contenues dans les formules. La seule intégrale, / "Trdt, four- 


5 
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t 
D D , . Ñ 10 A 
nira immédiatement les variables w, p, a, €, x, v Ag É dt donnera à 


(t) 
son tour 9, d’où a, et en introduisant ces valeurs dans l'équation 


p 
36) ;. . 4 AA AA 
(26) r 1 — ecos(b — a)” 


oú l'on donnera aux quantités p, e el « leurs valeurs au temps déterminé /, 
cette équation sera celle d'une conique qui représentera la trajectoire réelle 
d'une maniére plus exacte que la conique de la premiére approximation. 

Les données approchées de la trajectoire fourniront de méme immédiate- 
ment à l'aide de P, les trois intégrales contenues dans », 4, y, et, par consé- 
quent , les valeurs de ces variables. 

La seconde approximation de la trajectoire étant ainsi connue, on s'en 
servira de la méme maniére pour en déterminer une troisiéme, et ainsi de 
suite , cinq intégrales suffisant dans tous les cas. 

Si les forces F, T, P émanent d'autres points, décrivant également des 
trajectoires, il sera nécessaire de connaitre aussi pour la détermination de ces 
forces, les trajectoires approchées de ces autres points. Les formules serviront 
A trouver ainsi, avec une exactitude de plus en plus grande, les trajec- 
toires décrites par des points en nombre quelconque soumis à leurs actions 
mutuelles, que ces trajectoires soient de méme nature ou de natures diffé- 
rentes, et quelles que soient d'ailleurs leurs situations respectives. 


$ 22. Pour donner un exemple de l'application de nos formules à une 
trajectoire différente d'une conique, supposons (fig. 5) une masse m’ se 
mouvant sous l'action d'une vitesse initiale V et sollicitée par une force 
constante T, trés faible, perpendiculaire au rayon vecteur 7 émanant d'un 
centre O, et située dans le plan de V et de r (plan de l'orbite) (**). 


* F 2 pen , 
(*) La formule DU) a == donne déjà une valeur approchée du grand axe sans 
2Fr + v? 


aucune intégration, si l'on se contente de la valeur déjà connue de v. 
(**) Ce cas rappelle celui d'un póle d'aimant se mouvant sous l'action d'un courant vertical 
trés faible, o. 
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La première approximation de la trajectoire sera donnée par la ligne droite 
méme, CE, que décrirait m’, si T n'existait pas. Soit AOB une ligne fixe, 
origine des angles œ du rayon vecteur dans le plan de l'orbite; 0C = / une 


NS 
perpendiculaire à CE, 8 l'angle AOC et z l'angle de r et de V. 


Fic. 5. 


L'équation de la ligne droite CE sera, en supposant m! en C pour t = 4,, 


r=V B + VA — ty, 
l 


MNT eher? 


les paramètres étant /, 8, V, t,. 
En fonction de, ®, r, z, V, t on trouve facilement 


l —rsins, 


r coso 


3 


t — t = + 


tg(^-— £) = + coto, 
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équations qui serviront à déterminer /, £, 8 quand on connaitra V, =, r, d 
en fonction de £. 

Or, les bur générales (17), (18), (19) donnent dans le cas actuel, 
pour les valeurs a‘, V’, éi relatives à la nouvelle approximation de la trajec- 
toire : 


VI+T/ rdt Ve TK rdt 
一 — O ient Dees 
(ku dues lg c dr AC CN V que 
r — 
dt 
el 
t dt 
LL 7 E n ET v=p S LW T rd. 
t ti 


A l'aide de la première approximation, on a 


E dr an. 
r=VB + V(t LB. qe CR 


, — 1 
S rt f. VT + V — uy - dt = Pd Oy 


en posant 
| CUP RUNE, 
Eé og oe mp et 


Pour / — 4, n = 1. Donc 


On en déduira pour les nouvelles valeurs z' et Y). 


p 
Vi + T 
POS + xl Lh d 
EE Pa —— 
Pre y° 
ETS 
RTE 


sins VE + V*(t — 1} 
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A l'aide de ces valeurs et des équations 


l =r sino’ 


L iW T 
t— Ae 


3 
tg(9' — f') = coto’, 


on calculera les paramètres variables /', f, et &' qui donneront pour un 
temps £ donné, la seconde approximation de la trajectoire réelle. Cette seconde 
approximation est la droite 
i 
r 一 —————— 
cos(b— f”) 


représentée (fig. 5) par m’, m,; et Pon obtiendra aisément soit analytique- 
ment, soit géométriquement comme le montre la figure, la position plus exacte 
m, de m', pour un temps déterminé quelconque. 

Pour éviter la résolution de l'intégrale double contenue dans la valeur (1 9”) 
de 9, résolution qui n'aurait ici aucun intérêt, on peut se contenter, pour 
obtenir &', d'introduire à la place de % dans l'équation 


tg( — $^) = coto”, 
sa valeur donnée par la premiére approximation. 


Ce qui précéde suffit pour montrer l'extension que peuvent prendre nos 
formules et la généralité de leur application. 


CAS PARTICULIER DE LA MÉCANIQUE CÉLESTE. 


8 23. Pour rendre les formules précédentes applicables au cas de la méca- 
nique céleste, et retrouver les formules de Wronski, il faut choisir pour le 
- point origine ou pôle du mouvement, la masse centrale m du système; T et P 
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seront les résultantes des composantes des forces perturbatrices, par exemple 
des attractions de masses étrangères m, m!!!, ..., perpendiculaires au rayon 
vecteur de l'orbite de m’ dans le plan de cette orbite, et perpendiculaires à ce 
plan, dans le mouvement relatif de m’ par rapport à m. Quant à F, elle 
deviendra 

m + ml 


M 
+ R = — — + R (en posant M = m + m’), 
r r 


R étant la résultante des composantes radiales des forces perturbatrices (par 
exemple, la somme des attractions de »/, m" ... suivant le rayon vecteur) 
dans le mouvement relatif; il suffira donc de remplacer F par 一 去 十 R, pour 
obtenir les formules de la mécanique céleste. 

Nous allons réunir ces formules dans le tableau suivant : 


t 
(13) wp = (who Trd, Dë o=o f Da 


———— nÁÁ—sÁMX 
M —Rr? 
(44) w= ——————. 
(wp) +) Trdt 
t 2 
Iren, + Së pl 
(15) p= E a EEN 
M — Rr? y 
t 
ep, +S Trdt y 
17) tgo = LIRE Penn clita 
Ke r dr 
t t 
` (rbl +. 7. Trdt (wp) s d Trdi 
8 — 一 一 一 a] == UN T I OMBRE... o 
CED 1 r sinc : (46) (v) r 2 
T nw? M — Ry? 
(22) =——— où (4) = LÀ (21^) a w 
, 9; — rU D 7 
P — — 9Rr — v* 
r sinto r 
at |, d d 
(28) ¡e = Vi bé (29) e= (e), si 上 eos 一 ajo + esin?(b — a) 75 
a? " ' p WN 


www.rcin.org.pl 


DE LA MÉTHODE DE WRONSKI. 59 


d-(pw) 1 dp 


(25) tg(d — x) = i 
pw e p 


— (50) a (ain +, / sin —afeos(o-a) 
igo " 
(26) r — 4 


1 — e cos (D — a) 


Pdt 


v sins 


t 
(53) = (+ eost — + 9) de, (52) de = 


dud uat eu de 
(36) del $. f A do. 


sin Y 
(4) 


t f 4 
(57) A sin(x— 4-9)tg- de. 


Ce qui a été dit sur l'usage de ces formules dans le cas général s'applique 
identiquement et sans difficulté au cas actuel. 


§ 24. Dans l'établissement des formules précédentes, nous avons considéré 
comme fixe le plan auquel, par les coordonnées 4, n, est rapporté le plan 
variable de l'orbite. Il peut cependant se présenter des cas où c'est à un 
plan variable lui-méme que ce plan de l'orbite est rapporté; c'est ce qui 
arrive, par exemple, dans la théorie du mouvement de la lune, op l'on prend 
pour plan origine l'écliptique, c'est-à-dire le plan variable de l'orbite terrestre. 
11 faut done connaitre les coordonnées angulaires du plan variable et de la 
ligne fixe de l'orbite par rapport au plan fixe, ou plan invariable du systéme 
tout entier, — coordonnées que nous désignerons par 


H, V, X, — 


en fonction de coordonnées de cette orbite variable par rapport au plan 
origine variable lui-même, ou plan ¿ntermédiaire, coordonnées désignées 
par 


eps 95,9 
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et des coordonnées du plan intermédiaire par rapport au plan fixe, coordon- 
nées désignées par 
(ai (+) (x). 
Dans la figure ci-après (fig. 6) : 
mFB est le plan fixe, mBA le plan intermédiaire, 
mF sa ligne fixe; mA sa ligne fixe; 


mCD le plan de Porbite, 


mE sa ligne fixe. 
Il s'agit d'obtenir : 


H, v, X en fonction de y, y, x et de (7), (4), (x). 


Fic. 6. 


On aura, d'abord, dans le triangle DCB , 


cos H = cosy cos (4) — sin» sin (y) cos DB 
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et, comme 
cos H = cos y cos (1) — sin y sin (y) cos | (x) 一 (4) + ¿7 


(38) AE da 
sin H — VW 1 — cos H, 


le signe 十 pour le nœud ascendant. 


Pour 
Y =C0F, 
on a 
cos Y — cos BF cos BC — sin BF sin BC, 
sin Y — sin BF cos BG + cos BF sin BC 
où : 
BF — (y), 
sin BC sinDB ` sim[ (o) — (4) + y] 
sin y da sinH sin H 


dans le triangle DBC. 
Dans ce méme triangle, on a 


cos BC sin H — eos [ (x) — (y) + y lain cos( 人 + sin (y) cos». 
Par conséquent 


zem 


cos Y 一 


— (y) + 4] siny cos(») + sin (y) cos y 


L^ O 


sin Y 一 | cos [(x) — (4) + p ] sin y eos (y) + sin (y) cos» | 
| e z HTW (y) + 4] sin y. 
ou bien 
cos Y = len eos ( +) jeos[ (x )—Q) ++] cos (7) + sin (4) cot sl — sin Di sin | (x) Adel 
(59) 


sin y 


sin Y = SS? sin (4) eos| (x) — (y )-2] cos (y) + sin (x) cot RS + Cos (y) sin (69—( ed). 


6 
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Enfin, on aura pour 
X — CF + EC 
sin X = sin v cos(X — v) + cos Y sin(X — v), 


(40) x : 
cos X = cos y cos (X — V) — sin Y sin(X — v), 


expressions dans lesquelles 


sin (X — Y) = sin (x, — y) cos DC + cos(x — y) sin DC, 
cos (X — Y) = cos(x — y) cos DC — sin(x — y) sin DC. 


Mais, dans le triangle DCB, on a 


sin DC — sin) sin [69 — (9 + +] 
sin H 
cos DC — cos | (2) — (4) + +] duid cos y + SIN y eos 9) 
sin H 
Par conséquent, 
NS sin(#) 
sin (X—Y)— m jsin( x— 4) (sin 1 COU(#) + cos| (9) — LOC cles s | + eos(,— 0 sin| (x) rdl, 


ceos(X—1) = Se: 


sin(y) | cos(x, — y) | sin y cot(y) +cos| (x)— (4)-+ 4 ]eosy | —sin(x—)sin[ (7)—(2) 4- 9] : 
On obtiendra ensuite sin X el cos X par substitution dans les expressions 
données plus haut (7). 


$ 25. Si l'on connait, par exemple dans la théorie de la lune, les coor- 
données (4), (n), (x) de l'écliptique par rapport au plan fixe et les coordon- 
nées y, n, x de l'orbite lunaire par rapport à l'écliptique, les équations précé- 
dentes feront connaitre À coordonnées Y, H, X de l'orbite lunaire par 
rapport au plan fixe. (4), (n), (x) sont données par des équations semblables 
aux équations (35), SC (31) en fonction immédiate de la composante P' 
normale au plan de l'écliptique. 


(*) Tous ces calculs vérifient les formules (58), (59), (40) données sans démonstration par 
Wronski. (Réforme, ete., p. Lt. 
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p, n, x élant fonctions de ( $), (2), (x), Y, H, X, on obtiendra les varia- 
tions de d. n, x qui PERM jé variations (4), luis ) du plan de l'éclip- 
tique, en faisant varier seulement (4), (n), (x) ion K relations (38), 
(39), (40). On trouverait ainsi pour ces variations : 


Se Nx) YHX| 
à = fa | Ay) #) Xx) YHX? 
de = fs | AA) pibe 


fo fa, f; étant les trois fonctions résultant de la résolution des équations 
variées. Il ne faut pas oublier que, dans ces équations, V, H, X varient avec 
le temps, et ont des valeurs différentes pour chaque temps donné. L'intégra- 
tion des variations A, à, dx, de (1) à t, donnera les quantités Ai, »!, y’ dont 
y, i i ont varié pendant le temps / — (1) par l'influence des variations de 
(4), (x) , (x). D'ailleurs, les quantités y, », x, (35), (36), (37) obtenues pré- 
dnos sont, à une constante prés, celles dont les nda de l'orbite 
lunaire ont varié en supposant (4), (n), (x) constants. Il en résulte que les 
valeurs des coordonnées du plan de l'orbite lunaire par rapport à Pécliptique, 
seront égales à 

pci My 

YH Y, 


xx. 


On voit, cependant, que les caleuls nécessaires à la détermination précé- 
dente seraient trés compliqués et qu'ils exigeraient la considération des quan- 
tités auxiliaires (4), (1), (x), Y, H, X. 


8 26. On peut éviter l'introduction de ces quantités, et résoudre le 
probléme d'une maniére beaucoup plus élégante, en caleulant directement, 
pour l'écliptique, ses coordonnées par rapport au plan de l'orbite lunaire, pris 
pour terme de comparaison. Supposons que la terre soit prise pour point 
origine. L'écliptique sera le plan variable de l'orbite décrite par le soleil 
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autour de la terre. L'angle infiniment petit dy”, dont l'écliptique tourne autour 
du rayon vecteur du soleil dans le temps dt, sera donné par la formule 


P'dt 


DOULEUR TE THAI ES 
(41) d v'sins" 


P! étant la composante normale à l'écliptique des forces accélératrices qui 
sollicitent le soleil, v' sa vitesse réelle et w' l'angle de son rayon vecteur et de 
la tangente à Porbite qu'il décrit. Seient (fig. 7) 


m'S2, l'orbite lunairc; 


m" $2, l'orbite solaire. 


Nous supposons les forces P et P’ normales à leurs plans et dirigées dans 
| le sens de l'are ab qui mesure 
Pinclinaison mutuelle y. Si, 
par l’action de P”, Pécliptique 
se déplace de l'angle p', » 
deviendra »', y = QA de- 
viendra 4'— 6/A! en ayant 
toujours i 


Fic. 7. 


m"A=0, m'"A'--4', 


®' étant l'angle dont le soleil 
HUT s'est déplacé à partir de 
la ligne fixe mA, ou mA'. 
Enfin, y= m/6 + QA devient y! = mA! + Q'A’. 

Calculons dz, d, dy! pour ds n, Y! = d xy! — x. 


1. On a, en abaissant m’ perpendiculaire sur m/s, 


i EE Ee rie 
tgm"l=sini$3' tg4' et — —— — — ——, 
s 


inp sin y' 


d'ou l'on déduit, comme à la page 31, 


dy = — cot (53 siny cos y d, IG} 
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el, à cause de | 
d 6262' — + d-lQ', 


dy = — coti$2 sin m" Q cosy: del, 


d'oú encore (p. 31), 
(42) . . . . . . . dy — eosm"Q2dp' — — eos (9^ — de’. 
2. Pour caleuler dy, on a les relations 


d =m Q + y — DIR + y’, 


d’où 
dy = — d.m” Q. 
Or, 
sin m" Q’ sin m” Q 
SIN y sia Y ; 
d'ou 
sin sin (和 /一 
cos m" 2'd -m"62' 一 一 pe cosy’: dy”. 
sin? y 
Pour 
‘=, dy = dy 
et l’on a : 
sin (6' — y 
cos (5 — y)d- m” QR 一 mu c cos ($ — y) dp' 
ou 
sin (9' 一 
d "C2 ( H de”. 
ig # 
Donc : 
sin (和 "一 
(45) VUA. c duin CN 
1g y 


9. Enfin, pour dy, on a 


x — m'Q' + y — (x — y) + £20 + y”, 
dx' =d QY + dy’. 
Or, 
sin QG  sin(o'— y) 
"sinp. | aing 
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d’où 
d. 6262' aip sin (%'— y) 
Ur. CS ENSE 
Donc 
in(6'— in (9' — 
gy UA y, A y 
siny tg 4 
全 
= sin(' — y) ar Ar GEN 
sin 
et enfin 
CL M des DOCS E dx= + nE Y) 185 de" 


En intégrant maintenant, de (1) à t, les expressions (42), (43), (44), on 
obtiendra les quantités »”, 4, y! dont les coordonnées 5, 4, y ont varié par 
le fait des variations de position de Pécliptique. Ainsi 7 +7", 4 + V; x + y! 
seront, au temps £, les coordonnées de l'orbite lunaire par rapport à Péclip- 
tique. Ces coordonnées seront, en appelant (,),, Les Le leurs valeurs 
au temps E : 


m = (ma + SA cos(% — y + 0) dp 一 m cos (9^ — 4) de”, 


OI 


sin (x — y + 9) t sin(6' — y) 
AB) ui = de — ——— — — (lp, 
(45) du = (hi eie at p 4 ig p 


WM 


ya = (xa) + f. sin (y — y + 9) 187 de A sin(o' — y) tg de' Wi 
(1) (e) 

Il suffira donc, ici comme précédemment, de connaitre d'une facon 
approchée les coordonnées n, y, x de Porbite lunaire par rapport à Pécliptique 
fixe, l'angle éi et la variation absolue dp”, qui seront fournis par les formules 
connues, pour obtenir, par des substitutions successives, les valeurs de plus 
en plus exactes de », 4, x, en fonction du temps. 


$ 27. Nous possédons maintenant toutes les formules nécessaires pour 
déterminer par les méthodes algorithmiques connues, à l'aide d'approxima- 


(*) Données sans démonstration par Wronski. (Réforme, ete., p. c.xxxi.) 
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tions successives, toutes les circonstances du mouvement d'un globe, chaque 
fois que sa trajectoire réelle sera déjà connue d'une maniére approchée. 

Ce qui frappe dés l'abord dans ces formules, c'est leur absolue généralité. 
Que les trois forces R, T, P soient les composantes de l'attraction de globes 
extérieurs, ou dépendent seulement de la forme des masses attirantes, ou 
soient dues à la résistance d'un milieu, ou à une force répulsive de surface, ou 
aux déformations des globes, etc., ces formules donnent l'expression, sous 
forme d'intégrations finies, des paramétres variables de la conique, para- 
métres qui, substitués dans l'équation (26) de la trajectoire, donneront 
de celle-ci une représentation de plus en plus rigoureuse. 


$ 28. Comme M. Y. Villarceau l'a rappelé dans sa note sur les méthodes 
de Wronski, ce géométre opposait sa méthode, qu'il nomme méthode de 
l'ordre, à la méthode actuelle où Pon considère les forces exercées par les 


globes extérieurs comme des forces perturbatrices, méthode qu'il appelait 


méthode du désordre. 

Si, d'un cóté, l'on considére le mouvement conique d'un globe comme le 
mouvement qu'il devrait posséder, il est clair que toutes les forces qui trans- 
forment ce mouvement méritent le nom de forces perturbatrices. Si, d'un 
autre cóté, l'on regarde l'ensemble de toutes les forces motrices, comme assu- 
rant l'équilibre du systéme en l'assujettissant à des conditions périodiques, il 
n'y a aucune raison pour donner aux unes le nom de perturbatrices et aux 
autres celui de conservatrices. 

Toute la question serait néanmoins de savoir si l'ensemble de toutes ces 
forces concourt réellement à maintenir les conditions périodiques actuelles du 
système, ce qui est loin d’être prouvé, surtout en tenant compte des défor- 
mations des globes, qui doivent étre comprises aussi dans les formules géné- 
rales de Wronski. Il faut avouer, cependant, que le nom de forces perturba- 
trices est plus commode pour le langage qu'il n'est heureux au point de vue 
philosophique. 


8 29. Wronski a donné dans la Aéforme des mathématiques les expres- 
sions des paramétres variables dans la méthode du désordre, c'est-à-dire en 
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considérant les trois composantes générales R, T, P comme des forces 
perturbatrices. Nous les reproduisons ici pour deux raisons : 

4° D'abord, parce que la comparaison de cette méthode avec la précédente 
fera toucher du doigt l'esprit de cette derniére; 

2» Pour épargner à ceux qui voudraient les vérifier dans l'ouvrage et 
d’après la méthode de l'auteur, des calculs trés embarrassants et la substitu- 
tion de quelques raisonnements. 


8 30. Procédons à cette exposition. 

m' décrivant une conique autour de m sous la seule action réciproque de 
ces deux masses, on peut supposer, au temps /, les forces perturbatrices 
R, T, P appliquées à mi. chercher les variations qui résultent de ces forces 
seules, pendant le temps dí , dans le rayon vecteur », la vitesse v, l'angle z 
du rayon vecteur et de la tangente, l'angle dont le plan variable de l'orbité 
tourne autour du rayon vecteur, et ensuite en déduire, comme précédemment, 
par les propriétés des coniques, les variations des constantes de la trajectoire. 
C'est proprement le point de vue dans lequel les forces R, T, P sont pertur- 
batrices du mouvement conique primitivement établi. 

Des trois forces R, T, P, les deux premiéres seules, agissant dans le plan 
de l'orbite, exercent une influence directe sur la position de cette orbite dans 
son plan et sur ses dimensions, c'est-à-dire sur les trois paramètres 


Occupons-nous d'abord de leurs variations, et, pour cela, remarquons 
qu'ils sont fonction de trois autres quantités qui déterminent également d'une 
maniére compléte la trajectoire conique (puisque, en effet, la masse centrale 
in + m'est ici donnée), savoir : 


le rayon vecteur r correspondant à l'angle donné o, 
la vitesse totale v 


et l'angle w de la vitesse et du rayon vecteur. 


Or, les variations de ces trois derniéres quantités sont des fonctions immé- 
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diales des forces R, T qui étant, en effet, appliquées à la masse m', ont pour 
effet direct de faire varier sa vitesse, sa distance à m et sa direction. Il faut 
donc : 

1° A l'aide des relations connues entre <. 4, e etr, v, w dans le mouve- 
ment conique, exprimer les variations dz, da, de de o, a, e correspondantes à 
des variations arbitraires dr, dv, da de r, v, w; 

2» Chercher les différentielles dr, dv, dz de r, v, «+ produites par l'action 
de R, T pendant le temps dr | 

3° Poser ensuite dr — dr, dv = dv, da — dw et remplacer dans les expres- 
sions de ox, du, de. Ces variations seront alors les différentielles da, da, de 
produites pendant le temps dí par R et T et l'on n'aura plus qu'à intégrer 
pour obtenir «, a, e en fonction de /. 


$ 31. Les relations qui lient a, a, e à r, v, w ont été données précédem- 
ment. Ce sont les relations 


p 
SS)... Tm————., 
6) : 1 — ecosy 
(6) . . . (v)=w(1 —ecos y). 
(7) . . . v= wV1 + eê — 2e cosp 


auxquelles il faut joindre les relations également établies ou connues : 


(MT OL ee 
(10) .. (sf w=, p —a(1— e). 

io). v MO) == uestis EK? EE 
LUE id 


On tire de (10), 
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de (6) et (46), 


dr a 
(ve) 一 Wess = — we sing tgo, 


et les équations (8), (6), (7) deviennent, par substitution, 


(47) . . 7 1 — ecosg) = ol — e), 
(48) . . . e(cosp —singtgu)= 1, 
(49 . . . va(l —e)=M(1 + e? — 2e cos q), 


où n'entrent plus, outre l'angle p = ® — a, que les quantités a, e, a, r, v, w. 
(La quantité < est comprise dans cet angle < que nous avons conservé pour 
simplifier l'écriture.) | 


Remarque au sujet de la quantité h. Langle @ qui fixe la position du 
rayon vecteur 7, varie pendant le temps dí, et à cette variation dp de o 
répondent des variations (2r), (dv), (2) de r, v, «. 

Il en résulte que les variations totales et réelles de r, v, w sont égales à 


(dr) + dr, (du) + du, (ds) + du, 


en représentant par ð, à, d les variations dues à des forces étrangères. 
Si maintenant, l'on déduit des équations précitées les variations de a, e, a, 
on les obtiendra sous la forme 


da = A [ (3r) + dr] + A' [ (9v) + dv] + A" [(95) + ds], 
de — B [ (3r) + dr] + B' [ (9v) + 2v] + B" [ (25) + 95], 
da = C [ (3r) + dr] + C [(90) + w] + C" [(95) + ds], 


A, A', A”, B, B', B, C, C', C", étant des coefficients fonctions de a, e, a, 
r, V, w, 9. Or, en supposant nulle l'influence des forces étrangères, on 
aurait eu : 

o = A(dr) + A'(óv) + A"(óc), 

o = B(r) + B'(9v) + B" (2), 

o = C(dr) + C'(9v) + C" (25). 
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Par conséquent, 
da = Adr + A'óv + A" 0c, 


de — Bir + B'dv + B"óc, 
da = Cdr + C'óv + Cd. 


Il faut donc dans les équations (47), (48), (49) considérer (dr), (2v), (dx) 
comme nuls, ou, en d'autres termes, ® comme constant. Cette condition 
donne, en vertu de la relation ç = 9 — a, 


do — — da. 


8 32. Ceci étant bien établi, procédons d'aprés cette régle à la variation 
des équations (41), (48), (49). La variation de ces équations telles qu'elles 
sont données, fournit entre da, de, da trois équations du premier degré dont 
la résolution, et surtout la simplification ultérieure, donnent lieu à des calculs 
extrémement compliqués. On procéde plus simplement de la maniére suivante. 
On tire del'équation (47) : 
ail — e) 


ecosp=4— T 


En transportant cette valeur dans l'équation (49), elle devient : 


1—e 9 
v'a(1 -ej-u( EEN — dot UD —- wa — e) (2 —1) 
P 
ou bien, 
2a av? 
V — UM 


D'où, en variant, 
24a 2aàr vda + 2vadv 


r r? M : 
9 ^ 2a Iwa 
dal aiu e d a a 
de ce TM 
et, comme 
S ei 1 


r M a s 

2 Da? lo b : 

IUIS o e Ta m em dl + 二 óv. (H p / 
r A / 


/ 
/ 
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Pour obtenir & et de, varions maintenant les équations (47) et (48). Nous 
obtiendrons, pour (48) d'abord : 


ad 一 in edo 一 coso t Ud) 一 
e si co a =0 
e "T TS" costo < 
OU 

de | e sing 

— = e €08 o (t£ + tg aldo + de 

à 9 (tg? 5 ) de costo , 


où l’on a : 


y i : 4 -- ecoso e — COS? 
+ tgo = (89 一 —————— = ————— 
sli 8 ar e sine e sih ọ COS y 


1 1 + ecos? y — 2e cosy e + À — 2e cos; 
et — 一 1 + tgo — 1 + ——————À—————— — PR 


2 cin? 4 &in? 
| 608 万 e sine 6 an": 


el, par conséquent, 


de | e— cos 4 + e? — 2e cos 
(54) EE 
e sin e e sin 


La variation de (47) donnera 


r Ore ` t? cos g r'e sing 
dr 一 — da — — de + = 
a 


一 — De. 
q E ee 


D'où, à l'aide de l'équation (50), 


2avr , 2re j 1? cos? resiny 
D 一 — e + —— — — de — — h, c 
M 1-4 ett. éi at) ?” 


2a 
dr = — ór + 
r 


9wu2(1 一 2 
Zei SC 


2 
(52) . . . (r° coso — Qrae) de — rie sin ode = (1 — e?) a L — = dr M 
r ] 


cosp — e 


Multiplions l'équation (51) par er? sin, l'équation (52) par SE 
tons-les membre à membre, afin d'éliminer %. Nous obtiendrons : 


et ajou- 


K d SEH 
DIE 


r 


| (1 — ea (1 


sin ? 
Zo! — e)v cosp—e 


一 ————— x ——— . jy 


8 — 9rae) (e — 
|" PME: PE UE lae 
Mr sin y 


sin y 
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Le coefficient de 2r sera 


e— COS? aa e S Fo SEU 
r r — re cosp — 2ae coss + 2ae 
3 
a(1 — e?) E 1) 
e — COS ọ r 
RT T r(1 — ecos 9) — 2ae(cos o — e) 
2 
ali — éi ED 
e — coso r 
EE CES 
f ai =e) — due = À dae) 
e re 
2a 
a(1 — e) L 1) 
e — COS o T 
- - EE 
- vow 一 四 一 sa RN 2 
e 
2a 
a-a) 
€ — COS? Tu 
«t oh 3) 
. €— €089 
dics qu 


2(e — cos y) a4 — er av 
一 -一 一 一 .一 一 一 = e — eos y) == == 
Mr 20 24 
a(i — e? Eë den 
r r 
Mais, on a | 2a où 
| r M” 
done, ep coefficient se réduit à 
2(e — cos y) 
: 


Enfin, le coefficient de 2z sera 


(1 + e — re cosp)? sing rain: 4 + e 96 cos 


2a oT garg 2a 
[—e|—.— Eu EE leen 
ra( e 1) (1 le 1) 
et, à l'aide de (49), "sinp Call rain: 
e a aq a 
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On a donc , pour de, l'expression trés simple : 


/ 
| E v 
O T S TR ms LIED ue E di m 
a 


T v 


Quant à 35, l'équation (51) donnera 


sin 2 sin r sin? . À + e? — 9e cos 
UO MES UN IL. A A 
re ve ae(e — cos y) e(e — cos o) 


Le coefficient de d+ se simplifie de la manière suivante : 


r sino 1--e —2ecose  rsin*e—a — ae + 2ae cos 
ae(e — cos y) e(e — cos e) ae(e — cos e) 
r — a — r cos? yp — ae + 2ae cos o 


ae(e — cos) 
et, en remarquant que, d’après l'équation (47), r — a = re cosp — ae?, 


re cos p — 2ae* — r cos* o + 2ae cose 


ae(e — cos y) 


(r cos e — 2ae) (e— cos?) r cose — 2ae 


ae(e — cos y) ae 


de telle sorte que la variation da = — dp prend la forme définitive : 


sin < 2 sin y 2ae — r cos y W 
UM TE e ES cacher cei R. => dd: rr “ar” os (119 ^] 


Nous retrouvons ainsi dans les expressions (50), (53), (54) les varia- 
tions de a, e, «a telles que Wronski les donne sans démonstration sous le 
numéro (119) (p. cxxxv de la Réforme des mathématiques). 


$ 33. Il ne nous reste plus maintenant qu'à remplacer dans ces expres- 
sions, dr, dv, d par leurs valeurs en fonction de R et T. 


Variation de r. Le théoréme II nous donne immédiatement 


(68) 4, v 1.3 MS 4 A MITT 
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Variation de v. Sic est l'angle de la vitesse totale v et du rayon vecteur r 
au temps ¿, la résultante des nouvelles forces suivant la tangente sera 


R coso + T sinc, 


el la variation correspondante de » sera : 


(560. |... 5 boss, … dde — (B cosu.+ T sinc) dt. 


Variation de z. La décomposition de R et T suivant la perpendiculaire à 
la tangente, donne une résultante 


R sino — T coso 


qui, pendant le temps di, agit à l'extrémité du bras de levier infiniment petit 
vdt pour faire tourner la tangente autour de son point de contact et diminuer 
l'angle de — dz. L'élément de vitesse 


(R sino — T cosz) dt 


gagné pendant di, dans la direction de cette résultante, égale le chemin par- 
couru, en mouvement relatif, par son point d'application — vdt - dz, divisé 
par le temps dí. On a donc 


vdt. do 
R sino — T coso dt 三 一 .一 一 一 一 
(R sine — T cos») Ge 
d’où : 
1 
(57) e Eee arn ARES de = da = — ~- (R sins — T cosa) dt. 


$ 34. La substitution de (55), (56), (57) dans (50), (53), (54) donne 
alors, en effectuant immédiatement l'intégration de da — Au, da = ox, de — 9e, 
et en appelant 


(a), (e Lela 
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les valeurs de a, e, « pour une époque (1), 


"dav 
= (a), + 
(ole /5 M 


(0) 


Kau : (e) + Sr Are R coso + T sino M ura: R sinc — T coss) |dt, 
(58) (0 va 


\dt, 


2 sing Ces Zoe — r cos o d 
= Lë s c (R coso + T sinc) + S TRO C sins — T coso) |dt, 
( 


9 étant égal à D — a. 

Quant à la quantité 9 elle-même, où à l'angle que fait avec la ligne fixe du 
plan de l'orbite le rayon vecteur r, elle est égale au méme angle pris dans la 
conique variable à partir de la méme ligne fixe. Si p et w sont les paramétres 
variables de cette conique, on aura 


dé T. dt. 


Si donc (0), est la valeur de @ au temps (t), on aura 
t 
e w 
15: eget MUR de oc) f. P di, 
(2) 


ou encore, en fonction dea , e et M = m + m’, 
t —————————— 
A 1— e 
a 2 ares pecht ANE. A a 
D d 


puisque, dans cette seconde méthode, 


[MP Sp “一 V T. p=a(1 —e) et p» =V Mp =V Mall — e). 


$ 35. Les équations (47), (48), (49), (88), (59) renferment la solution 
du probléme par cette seconde méthode, du moins en ce qui concerne le 
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mouvement de m’ dans le plan de Porbite. Le mouvement de ce plan est pro- 
duit par la troisième force normale P; et les formules relatives à ce mouvement 
sont les mémes que dans la premiére méthode; ce qui devait étre, puisque, 
dans cette première méthode, on est également amené à considérer la force P 
seule, aprés avoir d'abord déterminé l'influence des forces R et T, et que cette 
force n'exerce aucune influence directe sur les éléments de l'orbite méme. Il 
n'y a donc pas lieu de les exposer de nouveau ici. 

Dans cette seconde méthode, comme dans la premiére, les forces R, T, P 
seront calculées en fonction du temps d'une maniére approximative, quand on 
connaitra par une premiére approximation la forme des trajectoires décrites 
par les centres attifants; il en sera de méme de la vitesse v, des angles ® et z 
et du rayon vecteur r. On introduira ces valeurs dans les intégrales et Pon 
obtiendra, en les résolvant, des valeurs approchées de a, e, « en fonction du 
temps /. A l'aide de ces valeurs, on en obtiendra de nouvelles plus exactes 
dev, w, 6, r, R, T, P..., dont on opérera de nouveau la substitution dans ` 
les intégrales, et ainsi de suite. Tel est du moins le procédé connu, qui, dés 
à présent, peut faire comprendre la praticabilité des formules. 


$ 36. On voit maintenant, d'une facon parfaitement claire, en quoi ` 
consiste la différence de la méthode de l'ordre de Wronski et de la méthode 
du désordre, où les forces R, T, P sont regardées comme des forces pertur- 
batrices. C'est que la conique variable est décrite, dans la seconde, sous 


l'influence des éléments 
Ds. Ye f, M, 


et dans la premiére, au contraire, sous l'influence des éléments, 


2. 
v, c, fr, —Fr?; 


C'est-à-dire que la masse constante M est remplacée, dans la méthode de 
Wronski, par une masse fictive variable — Fr°?. Ainsi, par exemple, pour la 
terre, dans la méthode des perturbations, l'ellipse variable sur laquelle elle 
ést supposée se mouvoir, est toujours décrite sous l'influence attractive de la 


8 
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masse du soleil; dans la méthode de Wronski, au contraire, l'ellipse variable 
de la terre est décrite sous l'influence d'une masse fictive, — Fr?, placée au 
centre du soleil, F étant la résultante, suivant le rayon vecteur r = ©, de 
toutes les forces qui agissent sur la terre, attractions du soleil et des pla- 
nétes, etc. 

Mais, quoique les orbites variables ne soient pas les mémes dans les deux 
cas, les coordonnées qu'elles fournissent pour la terre en chaque instant sont 
identiques, parce que ces coordonnées sont en chaque instant celles qui leur 
sont communes. Ces deux ellipses ont continuellement un plan commun (plan 
variable de l'orbite), une tangente commune, une vitesse v, un angle w, un 
rayon vecteur commun, et ces quantités communes sont également celles de 
l'orbite réelle de la terre au méme instant. 

Nous comprenons ainsi cette déclaration de Wronski, reproduite déjà par 
M. Villarceau, au sujet de la théorie lunaire : « Tout ce que nous pouvons 
» dire ici concernant la fausse théorie de la science actuelle, c'est que l'orbite 
» dela lune, qui y résulte des prétendues perturbations causées par le soleil , 
» n'est nullement identique avec l'orbite variable que découvre la vraie 
» théorie, par l'influence téléologique du soleil. » Nous comprenons aussi, 
comme l'a déjà remarqué M. Villarceau, que la concordance entre les deux 
systémes, annoncée par Wronski, ne doit s'entendre absolument que des 
résultats définitifs, c'est-à-dire de l'expression des coordonnées en fonction 
du temps. Les formules et l'esprit des méthodes sont entiérement différents, 
et, quant à la simplicité des caleuls, indépendamment du nombre des intégra- 
tions qui est de cinq dans la méthode de Wronski, et de 7 (58 (trois équa- 
tions), 59, 35, 36, 37) dans la deuxiéme méthode, il suffit de comparer des 
deux parts les expressions à intégrer, pour étre convaincu de l'avantage de la 
première. | | | 

La généralité et la variété des applications sont aussi tout en sa faveur. Ce 
point a été développé plus haut. La seconde méthode ne s'applique qu'à la 
circonstance particuliére oú la trajectoire réelle est à peu prés une conique. 


S 97. Enfin, la comparaison de ces deux coniques variables et de la tra- 


H 


jectoire réelle, nous conduit de plus à une conséquence trés importante, 
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contenue déjà dans le principe A; c'est que la conique de Wronski coincide 
mieux avec cette trajectoire réelle. 

En effet, si l'on représente par r, 7', r!! les rayons vecteurs de la trajec- 
toire réelle, de la conique de Wronski et de la conique actuelle, on aura, 
d'abord, en chaque instant, 


Mais cherchons les expressions des dérivées secondes 


dr’ dr” 


dë dë 


Si M est la somme des masses atlirante et attirée, et F la résultante de 
toutes les forces suivant le rayon vecteur, on aura respectivement, dans les 
deux coniques, 


d'r' | (Er?) dq T MEN 
d ` är "tu di Mx io pow e 


dt 
dr"! M E psi 
dt)” 
Or, on a 


dt dt dt r ro’ 


5, NI o" étant les angles de position du rayon vecteur dans les trois trajec- 
toires ; et dans la trajectoire réelle, on a 


| de di 
Donc, 
gr’ dr e dr" TM dr M dir 
E AAA 


Par conséquent, si Pon développe les rayons vecteurs en fonction du 
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temps, les trois premiers termes du développement seront communs à la 
conique de Wronski et à la vraie trajectoire, tandis qu'il n'y en a que deux 
de communs dans la méthode actuelle. 

On peut douter que Wronski ait apercu cet avantage de sa méthode, 
puisqu'il dit (p. civ de la Réforme des mathématiques), que l'influence des 
forces R et T, nulle dans z, se manifeste dans la dérivée a et dans les 
dérivées des ordres suivants. 

Quant aux coefficients différentiels de 5, 0”, 0”, on trouvera facilement 


do de , de" 
dt di d 
d'al dep” d^o 


Geller sedeo ret emit 
do” do oft T4 9 do ü 
dë dë dt dq T d 


par où l'on voit que les deux premiers termes des développements de 0, ', éi 
en fonction du temps sont communs aux trois trajectoires , le troisième iden- 
tique seulement dans les deux coniques; mais que dans le quatrième terme 


de AT. reparait l'influence de la force perturbatrice R, spéciale à la conique 
de la seconde méthode. 


8 38. Nous terminons ici l'exposé de la nouvelle conception de Wronski. 
Cette premiére partie contient l'établissement des formules fondamentales, 
nécessaires et suffisantes pour les problémes de la mécanique céleste. Dans 
chaque cas, il conviendra de remplacer dans ces formules les trois forces 
R, T, P par leurs expressions relatives à ce cas; et il est trés important de 
remarquer que, gràce aux méthodes algorithmiques générales d'intégration, il 
ne sera pas toujours nécessaire de développer ces fonctions perturbatrices en 
séries, on pourra les introduire tout entiéres dans les formules générales, qui 
deviendront ainsi applicables à des cas queleonques d'inclinaison, d'excentri- 
cité, etc., des orbites. Il y a cependant des cas oü cette introduction des expres- 
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sions générales de R, T, P serait inutile, et, dans ces cas, leur développement 
en séries sera avantageux. 

En appliquant les méthodes algorithmiques de Wronski pour le développe- 
ment des fonctions à l'aide d'autres fonctions génératrices données, nous avons 
retrouvé (*) les expressions données par lui de R, T, P (**) dans le cas d'une 
planéte sphérique troublée par d'autres planétes également sphériques, ce qui 
fournit la deuxiéme approximation des mouvements planétaires. On peut donc, 
dés à présent, avec la correction indiquée, employer ces formules. Mais, 
comme ce n'est là, en définitive , qu'un cas particulier de formules générales, 
et comme surtout il fallait faire usage de procédés analytiques non encore 
démontrés, nous avons réservé ce développement pour la seconde partie 
de ce travail, qui comprendra, outre la démonstration de la /o? supréme 
( Technie , première partie), l'étude des nouveaux procédés algorithmiques du 
grand géomètre. 


(*) Avec une trés légère modification, indiquée par Wronski lui-même, p. cxxxi. 
(*) Réforme des mathématiques, p. Lex, 
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NOTE I. 


La loi (1) Gdt = — wd? lie la force accélératrice G au mouvement angulaire T et à la 
quantité w, qui, dans le mouvement conique, est une constante égale à la moyenne des 
vitesses linéaires extrémes. D'aprés Wronski, cette loi ne serait pas seulement, comme 
nous l'avons compris, une loi particuliére au mouvement conique, et qui ne subsiste pour 
une loi G différente de la gravitation, qu'à la condition de w variable, mais, d'une maniére 
absolument générale, la loi fondamentale de la dynamique. 

Sans nous arréter longtemps aux considérations philosophiques par lesquelles il cherche 
à établir cette loi, nous dirons qu'elle exprime l'équilibre permanent entre la force accélé- 
ratrice centrale G, quelle qu'elle soit, et la force inerte de séparation qui provient de 
l'inertie de la matière, et qui constitue une résistance au changement de direction. Cette 
seconde force, S, est donc toujours dirigée dans un sens diamétralement opposé à celui de 
la force G qui tend à produire ce changement de direction. Son action élémentaire Sd: doit 
étre proportionnelle à la variation de direction de la perpendieulaire à G, et cette variation 
est évidemment l'angle d? dont la direction méme de G a varié pendant le temps dt. 
w étant le facteur de proportionnalité, on a done, d’après Wronski, Sdt — wdọ, et l'équi- 
libre permanent entre les valeurs absolues de G et S, équilibre sans lequel, d'aprés lui, le 
mouvement libre sur la trajectoire serait impossible, exige que l'on ait toujours la relation 
fondamentale G = — S ou (1) Gdt —wdy (*). 

S serait, dans ces idées, la véritable force centrifuge. Dans le cercle, la constante w est 
la vitesse totale v, et l'on a 


vd y? 
村 
di 1 
puisque, r étant le rayon, on a 
dy 
v='r_— 
di” 


relation connue. 

Comme nous l'avons montré, en admettant les idées actuelles de la dynamique, la 
relation (1) où w est constant, n'est possible que si G est la loi de Newton. Il en résulte 
que l'admission de la loi (1) comme loi fondamentale et générale à laquelle doivent satis- 


(*) Épitre à S. M. l'Empereur de Russie, pp. 28-35. Metz, in-4e; 1851. 
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faire tous les mouvements centraux, est le renversement de toute la mécanique moderne, 
Wronski n'a pas reeulé devant cette conséquence; — nous nous contenterons de soumettre 
au leeteur les remarques suivantes : 


I. 


Gdt — — wdy, dit Wronski, est la loi fondamentale de la dynamique; il admet pourtant 
, : . dv . . . T SEN 
l'exactitude de la loi actuelle G — 77, mais comme cas particulier de la première. Et voici 
comment il en opère la déduction : il cherche ce que devient l'expression de la force 

L , . d MM r e . a H 
accélératrice G — — w T qui détermine le mouvement conique, dans le cas où la conique 
décrite se réduit à une ligne droite, comme, par exemple, serait le mouvement d'une 
comète dans une ellipse dont le paramètre p serait nul, ou l’excentricité égale à l'unité. 


Comme la vitesse extrème u est alors nulle, l'expression générale 


Vu + 2w(w — u)(1 — coss) 
de la vitesse v, devient 
pes V 2*(1 — cosp) 
ét, en différentiant v?, 
vdv = u’ sing - de. 


Mais, v n'est ici autre chose que la vitesse c, = — w sing, parallèle à l'axe de la 
éonique. On a done ess at ei 


— ww Sine * dv = w° sino: dọ 
et, par conséquent, 


dv = — w:o. 
Mais 
Gdt = — wde. 
Done 
Gdt = dv, 
ou enfin, 
iud 
dt 


« Et» dit l'auteur (page 458 des Prolégomènes), « c’est là notoirement la loi fonda- 
» mentale de toute la dynamique. Ainsi, cette loi fondamentale n'est rien autre qu'un cas 
» partieulier de notre loi supréme (1) de la mécanique céleste; et il se trouve avéré que 
» cette loi primitive, comme nous l'avons annoncé en la faisant connaitre, est généralement 
» la loi supréme de toute la mécanique. » 

Or, ce raisonnement est un cercle vicieux, puisque les expressions des vitesses v, et vs, 
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dont la composition fournit la vitesse v (comme il le dit à la page 257), ne sont autres 
que les intégrales de Zi ; dt et dt qui expriment les composantes (multipliées par dt) de 
la force accélératrice 6, suivant les deux directions rectangulaires. On se sert donc, pour 
démontrer la relation G =F, de cette relation elle-mème. 


II. 


Si Gdt — — wdọ est la loi générale de la dynamique, comment cette loi, quand on y 
considère G comme une inconnue, assigne-t-elle pour expression de G une fonction déter- 
minée,*la loi de Newton? 


. L'auteur donne pour la courbe générale de l'orbite que parcourrait un astre conduit par 
une force centrale G quelconque, une équation de la forme 


(A) + . . . reose(C + f Gsing- dt) r sing (C' + f G cosgdt) = C", 


C, C', C" étant trois constantes (*). 

Dans cette équation, dit-il, le rapport de dọ à dt est absolument indéterminé, de telle 
sorte que, pour une méme force anns G, en pourra représenter un nombre indéfini 
de trajectoires. Si l'on y assigne ya di “le rapport = S donné par (1), on gie sur l'équation 
des sections coniques, — ce qui est, en effet, exact. Si l'on y donnait à E ? la valeur qui 
résulte de la loi des aires de Képler, savoir i, k étant une constante, on trouverait 
toutes les trajectoires à forces centrales, dans lesquelles cette loi se vérifie, etc. 

Il ne faut cependant pas regarder de bien prés l'équation (A) pour reconnaitre qu'elle 
n'est autre chose que le principe des aires lui-méme. 

€ -f6 sing: dt, C' + f G cos - dt sont les vitesses v,, v, parallèles à deux axes 
rectangulaires ; r cos®, r sing les bras de levier correspondants, et, par conséquent, 


DT C089 — var sing = C" 
ou 


d; 


en appelant (va) la vitesse perpendiculaire au rayon vecteur; et ce résultat est absolument 
indépendant de la loi qui lie dg à dt. 


(*) Épitre à S. M. l'Empereur de Russie, p. 39. 
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Ainsi l'auteur, par une inconcevable méprise, admet le prineipe des aires et ensuite une 
loi fondamentale, Gdt 一 — wdo, dont l'existence générale n'est possible que si ce prin- 
cipe des aires n'est pas généralement vrai. 

D'ailleurs, quelle que soit la relation (A) dont on admet la inin si Gdt — — wdp 
est la loi fondamentale de la dynamique, il wy a pas à choisir le rapport; T il est toujours, 
quel que soit G, donné par cette loi méme, et par conséquent, d’après les conditions que 
s'impose Wronski lui-méme, sa loi fondamentale (1) n'est possible que dans les sections 
coniques, ou pour une force G agissant en raison inverse du carré des distances. 


IV. 


Dés qu'on admet pour exacte la loi G—2, la généralité du principe des aires s'en 
déduit, comme cela est connu. ll y a done contradiction à accepter la loi et à repousser ce 
prineipe. 

Comme, dans l'esprit de Wronski, G =F est un cas particulier de (1) Gdt = — wdy, 
il s'ensuivrait que l'on pourrait déduire ‘des cas particulier, une conséquence non 
contenue dans le cas général. 


V. 


Il résulte de ee qui précède que l'affirmation de Wronski, suivant laquelle le principe 
des aires ne serait vrai que dans le cas de la loi Newtonienne, non-seulement est en opposi- 
tion avec les résultats les mieux établis de la science depuis Newton, mais, ce qui est bien 
plus grave, en flagrante contradiction avec les principes de Wronski lui-même. 

Il faut donc considérer comme mal fondées ‘toutes les déductions par lesquelles il 
calcule en partant de la loi (1) dans l'Épiître à S. M. l'Empereur de Russie, pages 41-49, 
les courbes décrites sous différentes lois d'action G en raison d'une puissance de la 
distance (*). 。 


VI. 


: TA d r e 
On peut, il est vrai, considérer G — — wt comme une condition à laquelle serait sou- 


mise la force G dans tous les problémes, mais alors, w doit étre considérée, ainsi que nous 
l'avons fait dans le cours de ce travail, que la force totale soit centrale ou non, comme 


(*) Il trouve, par exemple, pour G en raison inverse du cube de la distance, et en raison directe de la distance, 
respectivement une ellipse dont le centre est le centre des forces, et une trajectoire dont les anomalies, y, ne 
sauraient dépasser certaines limites; — tandis que la premiére de ces trajectoires est une spirale et la seconde une 
ellipse, dont le centre est le centre des forces, ainsi que Newton le démontre déjà par des considérations purement 
géométriques. (Principia, lib. I, pr. X.) 
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une variable, et cette variable doit se plier, et à la relatien G — — we, et aux lois géné- 
rales de la dynamique, telles que le principe des aires. 

Nous avouons avoir hésité d’abord à reconnaitre cette grave erreur, dont les consé- 
quences se retrouveront peut-être malheureusement dans plusieurs autres travaux capitaux 
de l'auteur (*). Nous ne croyons pas, cependant, qu'on puisse nier la justesse des réflexions 
précédentes. Si elles sont exactes, eomme nous en avons la conviction, elles auront l'avan- 
tage de mettre en garde les adeptes enthousiastes du Messianisme contre les affirmations 
absolues de l'illustre savant. La leeture des ouvrages de Wronski]nous fait vivre avec un 
génie extraordinaire, qui s'est attaqué à presque toutes les connaissanees humaines; mais 
l'impétuosité méme de ce génie doit nous mettre en garde contre lui; travaillons donc à 
étudier les clartés, et consolons-nous des obscurités en nous rappelant que nul, quelque 
grand qu'il soit, n'échappe au préeepte d'Horace : Non omnia possumus omnes. 


(*) Voyez à ce sujet : Science nautique des marées, in-4». Paris, 1855. 3 
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NOTE Il. 


Si l'on prend l'équation conique 
CEE AEL aR T o = r[1 —ecos(0 — «)] — p, 
cette équation, en y faisant varier convenablement toutes les quantités, pourra représenter 
la trajectoire réelle. 

On obtiendra ainsi l'équation différentielle 

o = dr[1 — e cos(? — a)] + r [e sin(® — a) (do — da) — cos (+ — )de] — dp 

(dd, da, de, dp étant les différentielles obtenues par les lois données de la variation des 
paramétres), ou bien 


d 
1) v s o — p — + r[esin(® — ajdo — cos (t — 0)de] — dp — re sin (& — a) da, 


et, en divisant par p, ` 


_ dr _esin(t a)d& — cos(& — ajde — dp echt dL 


r 1 — e cos(b — a) 


en dénotant par L la valeur du second membre. 
Si maintenant l'on pose ` 


ob sin (9 — [x]) de — cos (+ — [x]) de 
ma 1 — e cos (+ — [«]) 


= — ds, 


en appelant [x] une valeur constante de æ, on aura 


dr cke, sin ($ — [4]) de — cos(6 — [x]) de 
1 — ecos(® — [a]) 
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et, en intégrant, 
Ces 
1—e cos (d — [«])' 


C étant une constante d'intégration et e la base des logarithmes népériens. 
Si l'on suppose z — o, € deviendra le paramétre [p] d'un mouvement conique ayant [a] 
pour longitude de l'aphélie. On a done 


iar hv a wo Hr M. 


Cette équation a la méme forme que l'équation (180) de Wronski (Réforme des mathé- 
matiques, page cum), mais elle en diffère en réalité par la valeur de la fonction z. En effet, 
Wronski considère z comme indépendante de «, du moins d'une manière immédiate, et 
cela parce que, aprés avoir remplacé la différentielle dæ en fonction de dp et dw dans 
l'équation variée (B), il regarde l'intégration des fonetions en « comme comprise dans 
l'intégration des fonctions en p et w. Voiei comment il s'exprime : 

« Ainsi ayant éliminé dans l'équation différentielle générale (l'équation (B)), la diffé- 
» rentielle de de l'aphélie en la remplacant par des fonctions différentielles équivalentes, 
» formées par d'autres variables, nous pouvons maintenant prendre l'intégrale générale 
» de cette équation, en considérant comme constante la quantité «, parce que, de la 
» maniére dont cette élimination de d« vient d'étre faite, en enlevant cette différentielle 
» avec tout ce qu'elle tenait de l'équation primitive (l'équation o=r[ 1— e cos(o — all — p) 
» il sera tenu compte de ia variation de cette quantité dans tout son ensemble, et cela 
» par l'intégration des fonetions différentielles qui se trouvent ainsi équivalentes à cet 
» ensemble de la quantité a. » 

Mais ce raisonnement n'est pas rigoureux, puisque ce n'est pas seulement de, mais 
aussi «, qu'il faudrait remplacer en fonetion de ces autres variables, pour que l'intégration 
des fonctions différentielles de ces autres variables tint compte de la variation compléte 
de x. On peut mettre ceci nettement en évidence par la comparaison des deux expressions 
du rayon vecteur données par les équations (26) et (C). 


On a 

36 PRENNE AA 
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et 
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À — e cos(® — Lall 


Ces deux équations représentent également la trajectoire plane réelle, quand on y rem. 
place les paramètres variables par leurs expressions en fonction des forces accélératrices. 
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On en déduit 


[p  4—ecos(5 — a) ` 


2 Tos cag SON 7 ET, 

p 1 — e cos (b — 4) 
, ¿sin (0 — [4]) de — rein — [»])de 
1 — e cos (è — [e])] — 


, 


ou 
e sin(^ — [«]) d? — cos(» — [x]) de 


— L+ 
1 — e cos(® — Lel 


a 
EI 
| 


qui est l'expression que nous avons trouvée plus haut, tandis que Wronski donne à dz 
seulement la valeur 


d 
E + esin(o — «) da, 


Il importe de remarquer, d'ailleurs, que les formes 


m [p]: s Lakat 
Tad "5 t7 a —[qges(— a) 


ou toute autre du méme genre, où E, E”, ... peuvent toujours être aisément exprimées, 
n'ont aueune valeur théorique spéciale; elles peuvent seulement avoir une valeur pratique 
dans les cas où la détermination des intégrales = ou =" présenterait des facilités partieu- 
liéres. C'est ainsi que le rayon vecteur réel serait lié au rayon vecteur [r] d'une conique 
constante par la relation 


L'unique avantage de cette forme est de ramener immédiatement r à [r] quand la fonetion 
perturbatrice £' est nulle. 
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